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Notation 



We adopt Feynman's notation and denote a spacetime point by x = x^ = (t, x), fi = 
0, 1, • • • , d, with d the number of space dimensions, while the energy fco and momentum 
k of a particle will be denoted by k = k^ = (k , k). The time derivative <9 = d/dt 
and the gradient V are sometimes combined in a single vector = (do, — V). The 
tilde on is to alert the reader for the minus sign appearing in the spatial components 
of this vector. We define the scalar product k ■ x — k^x^ = k n t — k • x and use 
Einstein's summation convention. Because of the minus sign in the definition of the 
vector dfj, it follows that d^a^ = doao + V • a, with a p an arbitrary vector. 
Integrals over spacetime are denoted by 



When no integration limits are indicated, the integrals are assumed to run over all 
possible values of the integration variables. 

Natural units fi = c = fee = 1 are adopted throughout. 




x 



while those over energy and momentum by 




dfc d d k 
2^(27r) d ' 



©Amstex 



1 



Chapter 1 



Prelude 



Continuous quantum phase transitions have attracted considerable attention in this 
decade both from experimentalists as well as from theorists. (For reviews see Refs. 
[jj], ^J, ^ These transitions, taking place at the absolute zero of temperature, are 
dominated by quantum and not by thermal fluctuations as is the case in classical finite- 
temperature phase transitions. Whereas time plays no role in a classical phase transi- 
tion, being an equilibrium phenomenon, it becomes important in quantum phase tran- 
sitions. The dynamics is characterized by an additional critical exponent, the so-called 
dynamic exponent, which measures the asymmetry between the time and space dimen- 
sions. The natural language to describe these transitions is quantum field theory. In 
particular, the functional-integral approach, which can also be employed to describe 
classical phase transitions, turns out to be highly convenient. 

The subject is at the border of condensed matter and statistical physics. Typical 
systems being studied are superfiuid and superconducting films, quantum-Hall and re- 
lated two-dimensional electron systems, as well as quantum spin systems. Despite the 
diversity in physical content, the quantum critical behavior of these systems shows sur- 
prising similarities. It is fair to say that the present theoretical understanding of most 
of the experimental results is scant. 

The purpose of these Lectures is to provide the reader with a framework for study- 
ing quantum phase transitions. A central role is played by a repulsively interacting 
Bose gas at the absolute zero of temperature. The universality class defined by this 
paradigm is believed to be of relevance to most of the systems studied. Without im- 
purities and a Coulomb interaction, the quantum critical behavior of this system turns 
out to be surprisingly simple. However, these two ingredients are essential and have to 
be included. Very general hyperscaling arguments are powerful enough to determine 
the exact value of the dynamic exponent in the presence of impurities and a Coulomb 
interaction, but the other critical exponents become highly intractable. 

The emphasis in these Lectures will be on effective theories, giving a description 
of the system under study valid at low energy and small momentum. The rationale 
for this is the observation that the (quantum) critical behavior of continuous phase 
transitions is determined by such general features as the dimensionality of space, the 
symmetries involved, and the dimensionality of the order parameter. It does not depend 
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on the details of the underlying microscopic theory. In the process of deriving an 
effective theory starting from some microscopic model, irrelevant degrees of freedom 
are integrated out and only those relevant for the description of the phase transition are 
retained. Similarities in critical behavior in different systems can, accordingly, be more 
easily understood from the perspective of effective field theories. 

The ones discussed in these Lectures are so-called phase-only theories. The are the 
dynamical analogs of the familiar 0(2) nonlinear sigma model of classical statistical 
physics. As in that model, the focus will be on phase fluctuations of the order param- 
eter. The inclusion of fluctuations in the modulus of the order parameter is generally 
believed not to change the critical behavior. Indeed, there are convincing arguments 
that both the Landau-Ginzburg model with varying modulus and the nonlinear O(n) 
sigma model with fixed modulus belong to the same universality class. For technical 
reasons a direct comparison is not possible, the Landau-Ginzburg model usually being 
investigated in an expansion around four dimensions, and the nonlinear sigma model 
in one around two. 

In the case of a repulsively interacting Bose gas at the absolute zero of tempera- 
ture, the situation is particular simple as phase fluctuations are the only type of field 
fluctuations present. 

These Lectures cover exclusively lower-dimensional systems. The reason is that it 
will turn out that in three space dimensions and higher the quantum critical behavior is 
in general Gaussian and therefore not very interesting. 

Since time and how it compares to the space dimensions is an important aspect of 
quantum phase transitions, Galilei invariance will play an important role in the discus- 
sion. 
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Functional Integrals 



In these Lectures we shall adopt, unless stated otherwise, the functional-integral ap- 
proach to quantum field theory. To illustrate the use and power of functional integrals, 
let us consider one of the simplest models of classical statistical mechanics: the Ising 
model. It is remarkable that functional integrals can not only be used to describe quan- 
tum systems, governed by quantum fluctuations, but also classical systems, governed 
by thermal fluctuations. 



The Ising model provides an idealized description of an uniaxial ferromagnet. To be 
specific, let us assume that the spins of some lattice system can point only along one 
specific crystallographic axis. The magnetic properties of this system can then be mod- 
eled by a lattice with a spin variable s(x) attached to every site x taking the values 
s(x) = ±1. For definiteness we will assume a d-dimensional cubic lattice. The Hamil- 
tonian is given by 



Here, x = aa^e^, with a the lattice constant, Xi integers labeling the sites, and ej 
(i = 1, • • • , d) unit vectors spanning the lattice. The sums over x and y extend over the 
entire lattice, and J(x, y) is a symmetric matrix representing the interactions between 
the spins. If the matrix element J(x, y) is positive, the energy is minimized when the 
two spins at site x and y are parallel — they are said to have a ferromagnetic coupling. 
If, on the other hand, the matrix element is negative, anti-parallel spins are favored — 
the spins are said to have an anti-ferromagnetic coupling. 
The classical partition function Z of the Ising model reads 



with j3 = 1/T the inverse temperature. The sum is over all spin configurations {s(x)}, 
of which there are 2 N , with N denoting the number of lattice sites. To evaluate the 



2.1 Ising Model 




(2.1) 




(2.2) 
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partition function we linearize the exponent by introducing an auxiliary 0(x) at each 
site via a so-called Hubbard-Stratonovich transformation. Such a transformation gen- 
eralizes the Gaussian integral 



exp(-i/3J s 2 ) = ^-^^exp(-i/3J-V 2 + /?<^), 
where the integration variable <f> runs from — oo to oo. The generalization reads 



exp 



±3^J(x,y)s(x) S (y) 



(23) 



(2.4) 



J- X (x,y) 0(x) 0(y) + /?]T ^(x)s(x) 



Here, J _1 (x, y) is the inverse of the matrix J(x, y) and we ignored — as will be done 
throughout these notes — an irrelevant normalization factor in front of the product at 
the right-hand side. The equation should not be taken too literally. It is an identity 
only if J(x, y) is a symmetric positive definite matrix. This is not true for the Ising 
model since the diagonal matrix elements J(x, x) are all zero, implying that the sum 
of the eigenvalues is zero. We will nevertheless use this representation and regard it as 
a formal one. The partition function now reads 



z= j2 n/ ^wexp 



{s(x)} 



±/3^T j ~W) 0(x) 4>(y) + /3]T 0(x) s(x) 



(2.5) 

The spins are decoupled in this representation, so that the sum over the spin configura- 
tions is easily carried out with the result 

Z = Uf #(x)ex P Ul/^J- 1 (x,y)0(x)0(y)+^ln{cosh[/30(x)]}) , 

x \ x,y x / 

(2.6) 

ignoring again an irrelevant constant. 

The auxiliary field 0(x) is not devoid of physical relevance. To see this let us first 
consider its field equation: 



(x)=£j(x,y),s(y) 



(2.7) 



which follows from (2^). This shows that the auxiliary field </>(x) represents the effect 
of the other spins at site x. To make this more intuitive let us study the expectation value 
of the field. For simplicity, we take only nearest-neighbor interactions into account by 
setting 



J(x,y) 



J if site x and y are nearest neighbors 
otherwise, 
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with J positive, so that we have a ferromagnetic coupling between the spins. The 
model is now translational invariant and the expectation value (s(x)) is independent of 



x: 



( S (x)) = M. 



(2.8) 



We will r efer to M as the magnetization. Upon taking the expectation value of the field 
equation ( 2/7 ), 

<0(x)) = 2dJM, (2.9) 

where 2d is the number of nearest neighbors, we see that the expectation value of the 
auxiliary field represents the magnetization. 

A useful approximation often studied is the so-called mean-field approximation. It 
corresponds to approximating the integral over </>(x) in (2.6) by the saddle point — the 
value of the integrand for which the exponent is stationary. This is the case for (/>(x) 
satisfying the field equation 



^J- 1 (x,y)0(y)+tanhLS0(x)]=O. 



(2.10) 



We will denote the solution by m f. In this approximation, the auxiliary field is no 
longer a fluctuating field taking all possible real values, but a classical one having 
the value determined by the field equation ( 2.10 ). Being a nonfluctuating field, the 
expectation value (</> m f (x)) = m f (x), and ( 2.10| ) yields a self-consistent equation for 
the magnetization 

M = tanh(2d/3 JM), (2.11) 



where we assumed a uniform field solution and invoked Eq. (2.9). It is easily seen 
graphically that the equation has a nontrivial solution when 2d(3 J > 1. If, on the other 
hand, 2d(3J < 1 it has only a trivial solution. It follows that 



fa 1 = 2d J 



(2.12) 



is the critical temperature separating the ordered low-temperature state with a nonzero 
magnetization from the high-temperature disordered state where the magnetization is 
zero. 

Let us continue by expanding the Hamiltonian in powers of <\>. To this end we note 
that the term ln[cosh(/30)] in (2.6) has the Taylor expansion 



ln[cosh(/3</>)] 



12 



(2.13) 



Before considering the other term in (|2.6|), ^ x J _1 (x, y) 0(x) <fi(y), let us first study 
the related object ^ x J(x, y) s(x) s(y) which shows up in the original Ising Hamil- 
tonian (2.1). With our choice (pTq) of the interaction, the Taylor expansion of this 



object becomes 



J2 A*, y) *W s(y) = J S W ( 2d + ' 



; 2 V 2 



•) S (x) 



(2.14) 
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neglecting higher orders in derivatives. From this it follows that 



J-^x, y) 0(x) 0(y) = J- 1 £ 0(x) f i - ^a 2 



2 V 2 



and the partition function (2.6) becomes in the small-0 approximation 



z = \{j #w e ~ 



0H 



(x), (2.15) 



(2.16) 



with H the so-called Landau-Ginzburg Hamiltonian 



H 



E 



8d 2 J 



(V</>) 2 



-(— 

2 \2dJ 



12 y 



(2.17) 



The model has a classical phase transition when the coefficient of the 2 -term changes 
sign. This happens when = l/2dJ in accord with the conclusion obtained by in- 
specting the self-consistent equation for the magnetization ( {2.1 1[ ). 

In the mean-field approximation, the thermal fluctuations around the mean-field 
configuration are ignored, so that becomes a nonfluctuating field. The functional 
integral J\ x J d</>(x) is approximated by the saddle point. 

For future reference we go over to the continuum by letting a — > 0. To this end we 
replace the discrete sum J^i by tne integral a~ d f , and rescale the field <j>(x), 



(x)-0'(x) = 



0a 



2-d 



Ad 2 J 



(x), 



(2.18) 



such that the coefficient of the gradient term in the Hamiltonian takes the canonical 
form of i. In this way the Hamiltonian becomes 



0H = 



(V0) 2 



1 



X (j) 



(2.19) 



where we dropped the prime on the field; the parameter ro and the coupling constant 
Ao are given by 



'■{) 



^(0o-P), 



An = 



J 2 /3 4 



The partition function now reads 



Z = / D^e 



-0H 



(2.20) 



(2.21) 



where the functional integral J denotes the continuu m limi t of the product of inte- 
grals Yl x J di/)(x). The last two terms in the integrand of ( 2.19 ) constitute the potential 

V(0), 

' ' ' ' ' (2.22) 



1 



In Fig. Z. 1 , the potential V(<f>) is depicted in the high-temperature phase where ro > 0, 
and also in the low-temperature phase where ro < 0. The minimum of the potential 
in the low-temperature phase is obtained for a value 0^0, whereas in the high- 
temperature phase the minimum is always at cf> = 0. 
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Figure 2.1: The potential V((j>) of the Ising model in the high-temperature (tq > 0) and 
low-temperature (rp < 0) phase. 



2.2 Derivative Expansion 



We are interested in taking into account field fluctuations around the mean field <p m f, 
which is the solution of the field equation obtained from (2.19). To this end we set 
<fi = 0mf + <t>, an d expand the Hamiltonian around the mean field up to second order in 



(V0) 2 + (r + §A ^ 



lf)<t> 2 



(2.23) 



where H m f denotes the value of the Hamiltonian (2.19) for <j) — cf> m [. Because of the 



change of variables, the functional integral J D<j) changes to J Dtf>. Since we neglected 
higher-order terms, the functional integral is Gaussian and easily carried out. The 



partition function (|2.21|) becomes in this approximation 

"(V0) 2 + (r + iA o 2 nf )^ 



Dcf> exp 



e -/3H mf Det -l/2 (p 2 + r() + 1 Ao^), 



(2.24) 



with the derivative p = — iV. The determinant represents the first corrections to the 
mean-field expression exp(— f3H m { ) of the partition function due to fluctuations. Using 
the identity Det(A) = exp [Tr ln(A)], we can collect them in the effective Hamiltonian 



l3H cS = iTrln[p 2 + r + |A ^ f (x)], 



so that to this order 



3 -/3(ffmf + -Hoff) 



(2.25) 
(2.26) 



As indicated, the mean field </> m f (x) may be space dependent. 
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We next specify the meaning of the trace Tr appearing in (2.25). Explicitly, 

PH eS = lJjn{[p 2 + r Q + |A <& f (x)] *(x - y)| y=x } • (2.27) 



The delta function arises because the expression in parenthesis at the right-hand side 
of (2.24) is obtained as a functional derivative of the Hamiltonian (2.23), 



S 2 (3H 



[p 2 + r + iAo^ f (x)] 5(x-y)| 



y=x' 



(2.28) 



which gives a delta function. Since it is the unit operator in fun ction space, the delta 
function may be taken out of the logarithm and we can write for ( 2.27 ) 



PH cS = i J In [p 2 + r + |A ^ f (x)] <5(x - y)| y=x 

= \J J e- lk x In [p 2 + r + \ A ^ 2 lf (x)] e* kx . (2.29) 



In the last step, we used the integral representation of the delta function: 



J(x) 



(2.30) 



shifted the exponential function exp(— ik • y) to the left, which is justified because the 
derivative p doe s not operate on it, and, finally, set y equal to x. We thus see that the 
trace Tr in ( 2.29 ) stands for the trace over discrete indices as well as the integration over 
space and over momentum. The integral J k arises because the effective Hamiltonian 
calculated here is a one-loop result with k the loop momentum. 

The integrals in (2.29) cannot in general be evaluated in closed form because the 
logarithm contains momentum operators and space-dependent functions in a mixed 
order. To disentangle the integrals resort has to be taken to a derivative expansion 
in which the logarithm is expanded in a Taylor series. Each term contains powers of 
the momentum operator p which acts on every space-dependent function to its right. 
All these operators are shifted to the left by repeatedly applying the identity 



/(x)p fl (x) = (p + »V)/(x)s(x), 



(2.31) 



where /(x) and g(x) are arbitrary functions and the derivative V acts only on the next 
object to the right. One then integrates by parts, so that all the p's act to the left where 
only a factor exp(— ik • x) stands. Ignoring total derivatives and taking into account 
the minus signs that arise when integrating by parts, one sees that all occurrences of p 
(an operator) are replaced with k (an integration variable). The exponential function 
exp(ik ■ x) can at this stage be moved to the left where it is annihilated by the function 
exp(— ik • x). The momentum integration can now in principle be carried out and the 
effective Hamiltonian be cast in the form of an integral over a local density 7i c ff : 



HeS- 



(2.32) 
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This is in a nutshell how the derivative expansion works. 



Let us illustrate the method by applying it to (2.25). When we assume (j> m { to be a 
constant field <f>, the effective Hamiltonian (2.25) may be evaluated in closed form: 



/3V, 



off 



M 



r(i-rf/2) 

d{An) d / 2 M d ' 



M 



r + ±A O 2 , (2.33) 



where instead of an Hamiltonian we introduced a potential V c ff to indicate that we 
are working with a space-independent field cf>. To obtain the last equation, we first 
differentiated ln(fc 2 + M 2 ) with respect to M 2 and used the dimensional-regularized 
integral 

1 T(a-d/2) 1 



k (k 2 + M 2 ) a (47r) d / 2 r(a) ( M 2^- d / 2 



(2.34) 



to suppress irrelevant ultraviolet divergences, and finally integrated again with respect 
to M 2 . To illustrate the power of dimensional regularization, let us consider the case 
d = 3 in detail. Introducing a momentum cutoff, we find in the large-A limit 



1 



1 



87T^ 12-7T 



(2.35) 



where we ig nored irrelevant, ^-independent constants proportional to powers of A. We 
see that in ( 2.33 ) only the finite part emerges. That is, all terms that diverge with a 
strictly positive power of the momentum cutoff are suppressed in dimensional regular- 
ization. These contributions, which come from the ultraviolet re gion , cannot physically 
be very relevant because the simple Landau-Ginzburg model ( 2.19 ) stops being valid 
here and new theories are required. It is a virtue of dimensional regularization that 
these irrelevant divergences are suppressed. 

Expanded up to fourth order in </>, (2.33) becomes 



/3V C 



1 

12V 



3/2 



1 

Wir' 



-A n r, 



o'o 



1/212 



1 



\2 



1287ri/2 



(2.36) 



where the first term is an irrelevant ^-independent constant. These one-loop contribu- 
tions, when added to the mean-field potential 



/3V - \r^ 2 + ^A o 4 , 
lead to a renormalization of the bare parameters 



A = A °-i6^7^' r = ro 



8- 



Aor, 



1/2 



(2.37) 



(2.38) 



In the case </> m f is not a constant field, we write the mean field m f (x), solving the 
field equation, as (j) m f (x) = <f> + 0(x), where <f> is the constant field introduced above 



( J2.33[ ), and expand the logarithm at the right-hand side of ( |2.25| ) to second order in <fi: 

_J 2 1 

4"" P-+.U-'"" <S p2 + M 2 ^p 2 + M 2 ' 



,7/ "" = T^oTr 2 ^ M2 (2# + 2 ) ~ ^A^ 2 Tr 



(2.39) 
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with 



:= H c 



dV, 



cff 



(2.40) 



Moving the momentum operator p to the left by using ( 2.31 ), we obtain 

4 p z + M z 8 



p 2 + A/ 2 (p - iV) 2 + M 2 



(2.41) 

where we recall the definition of the derivative V as operating only on the first object 
to its right. Using the integral 



1 



1 



- 1 arctanf |q| 
k k 2 + M 2 (k + q) 2 + M 2 47r|q| \2M 



with q = — iV, we obtain for ( [2.4 1| ) 

(3H eS 



^A M(2# + <£ 2 ) - ^o<P ' 

1D7T OZ7T 



" 1 , 


f |q| V 


- — - arctan 






v2MA 



(2.42) 



(2.43) 



We note that only terms with an even number of derivatives appear in the expansion of 
this expression. The coefficient of the linear term is d(3V c s / d<fi, while that of the two 
quadratic terms independent of q is ^d 2 (3V c s/d4> 2 , as it should be. For Z we obtain 



z(4>) 



1 A 2 ^ 2 



(2.44) 



192tt M 3 ' 

Other terms involving higher powers of <f>, obtained from expanding the logarithm in 
(2.25) to higher orders, can be treated in a similar fashion. 
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Chapter 3 

Superfluidity 



A central role in these Lectures is played by an interacting Bose gas. In this chapter 
we wish to study some of its salient features, notably its ability to become superfluid 
below a critical temperature. We shall derive the zero-temperature effective theory of 
the superfluid state, and discuss the effect of the inclusion of impurities and of a l/|x|- 
Coulomb potential. Finally, vortices both at the absolute zero of temperature and at 
finite temperature are studied. 



3.1 Bogoliubov Theory 

The system of an interacting Bose gas is defined by the the Lagrangian [Q| 

C = <P* [ido - e(-<V) + Mo]<£ - AoH 4 , (3.1) 

where the complex scalar field 4> describes the atoms of mass to, e(— iV) = — V 2 /2m 
is the kinetic energy operator, and /i is the chemical potential. The last term with 
positive coupling constant, A > 0, represents a repulsive contact interaction. The 
(zero-temperature) grand-canonical partition function Z is obtained by integrating over 
all field configurations weighted with an exponential factor determined by the action 
S = f C: 

Z= I D0*D0e lS . (3.2) 



This is the quantum analog of Eq. (2.21 ) — the functional-integral representation of a 
classical partition function. 

The theory (3T) possesses a global U(l) symmetry under which 

<j>(x) -> e ia (j){x), (3.3) 

with a a constant transformation parameter. At zero temperature, this symmetry is 
spontaneously broken by a nontrivial ground state, and the system is in its superfluid 
phase. Most of the startling phenomena of a superfluid follow from this symmetry 
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Re(0) 



Im(<^>) 



Figure 3.1: Graphical representation of the potential (3.4) 



breakdown. The nontrivial groundstate can be easily seen by considering the shape of 
the potential 

V = - Alo |0| 2 + A o |0| 4 , (3.4) 



depicted in Fig. 3.1. It is seen to have a minimum away from the origin 
account for this, we shift by a (complex) constant cf> and write 



0. To 



(3.5) 



The phase field tp(x) represents the Goldstone mode accompanying the spontaneous 
breakdown of the global U(l) symmetry. At zero temperature, the constant value 



^ 2 = \y 

I Ao 



(3.6) 



minimizes the potential energy. It physically represents the number density of particles 
contained in the condensate for the total particle number density is given by 



nix) = \m\ 2 - 



(3.7) 



Because <f> is a constant, the condensate is a uniform, zero-momentum state. That is, 
the particles residing in the ground state are in the k = mode. We will be working in 
the Bogoliubov approximation which amounts to including only the quadratic terms in 
cp and ignoring the higher-order ones. These terms may be cast in the matrix form 



£(2) = Wm q (p,x)$, l> = 



(3.8) 
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with 

M (p,x)= (3.9) 

po-e(p)+Lio-U{x)-4\ \4>\ 2 -2\q4> 2 

~2\o4>* 2 -p -e(p)+iio-U(x)-4,\ \<p\- 

where U stands for the combination 

U(x) = d <p(x) + -L[V<^(x)] 2 . (3.10) 
2m 



In writing (3.9) we have omitted a term V tp containing two derivatives which is irrel- 
evant in the regime of low momentum in which we shall be interested. We also omitted 
a term of the form • j, where j is the Noether current associated with the global 
U( 1 ) symmetry, 

j = ^f Vf (3.11) 

This term, which after a partial integration becomes — ipW ■ j, is irrelevant too at low en- 
ergy and small momentum because in a first approximation the particle number density 
is constant, so that the classical current satisfies the condition 

V-j = 0. (3.12) 

The spectrum E(k) obtained from the matrix Mo with the field U set to zero is the 
famous single-particle Bogoliubov spectrum [Q], 

E(k) = + 2 w(k) 



/ e 2 (k)+4A o |0| 2 e(k). (3.13) 

The most notable feature of this spectrum is that it is gapless, behaving for small mo- 
mentum as 

£(k)~uo|k|, (3.14) 



with uq = y/ /j,o/m a velocity which is sometimes referred to as the microscopic sound 
velocity. It was first shown by Beliaev ^ that the gaplessness of the single-particle 
spectrum persists at the one-loop order. This was subsequently proven to hold to all 
orders in perturbation theory by Hugenholtz and Pines [g]. For large momentum, the 
Bogoliubov spectrum takes a form 

£(k)~e(k)+2A o |0| 2 (3-15) 
typical for a nonrelativistic particle with mass m moving in a medium. To highlight the 



condensate we have chosen here the second form in (3.13) where [1q is replaced with 

2A H 2 . 
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Figure 3.2: Graphical representation of the effective theory (3.22). The symbols are 
explained in the text. 



3.2 Effective Theory 

Since gapless modes in general require a justification for their existence, we expect 
the gaplessness of the single-particle spectrum to be a result of Goldstone's theorem. 
This is corroborated by the relativistic version of the theory. There, one finds two 
spectra, one corresponding to a massive Higgs particle which in the nonrelativistic 
limit becomes too heavy and decouples from the theory, and one corresponding to 
the Goldstone mode of the spontaneously broken global U(l) symmetry JlO[]. The 
latter reduces in the nonrelativistic limit to the Bogoliubov spectrum. Also, when the 
theory is coupled to an electromagnetic field, one finds that the single-particle spectrum 
acquires an energy gap. This is what one expects to happen with the spectrum of 
a Goldstone mode when the Higgs mechanism is operating. The equivalence of the 
single-particle excitation and the collective density fluctuation has been proven to all 
orders in perturbation by Gavoret and Nozieres JlT[] . 

Let us derive the effective theory governing the Goldstone mode at low energy and 
small momentum by integrating out the fluctuating field <t jl2[ ] . The effective theory 



is graphically represented by Fig. 3.2. A line with a shaded bubble inserted stands for i 
times the full Green function G and the black bubble denotes i times the full interaction 
T of the $-field with the field U which is denoted by a wiggly line. Both G and Y are 
2x2 matrices. The full interaction is obtained from the inverse Green function by 
differentiation with respect to the chemical potential, 

r-^. (3.16, 



This follows because U, as defined in ( 3.1C ), appears in the theory only in the combi 



nation /iq — U. To lowest order, the inverse propagator is given by the matrix Mq in 



(3.9) with U (x) set to zero. It follows that the vertex of the interaction between the $ 
and [/-fields is minus the unit matrix. Because in terms of the full Green function G, 
the particle number density reads 

n = l -tr [ G(k), (3.17) 



2 



we conclude that the first diagram in Fig. 



k 



3.2 



stands for —nil. The bar over n is to 



indicate that the particle number density obtained in this way is a constant, representing 
the density of the uniform system with U(x) set to zero. The second diagram without 
the wiggly lines denotes i times the (0 0)-component of the full polarization tensor, 
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IIoo, at zero energy transfer and low momentum q, 

ilimn O0 (0,q) = ~ lim tr / GT G (k , k + q). 

q^O 2 q->0 y fc 



(3.18) 



The factor 5 is a symmetry factor which arises because the two Bose lines are identical. 
We proceed by invoking an argument due to Gavoret and Nozieres [pT[ to relate the 
left-hand side of (3.18) to the sound velocity. By virtue of relation ( 3.16[ ) between the 
full Green function G and the full interaction V, the (0 0)-component of the polarization 
tensor can be cast in the form 



lim n 00 (0, q) 

q->0 



--limtr / G — — G(k , k + q) 
2 q^o J k dp 

i_d_ 

2 9/iq^u Jk 

On _ 1 ffl 



lim tr / G(fc ,k + q) 
Jk 



(3.19) 



where il is the thermodynamic potential and V the volume of the system. The right- 
hand side of (3.19) is n 2 n, with k the compressibility. Because it is related to the 
macroscopic sound velocity c via 

1 



(3.20) 

rrinc- 

we conclude that the (0 0)-component of the full polarization tensor satisfies the so- 
called compressibility sum rule of statistical physics [|TTj] 



lim n 00 (0, q) 

q^O 



n 2 K 



(3.21) 



Putting the pieces together, we infer that the diagrams in Fig. 3.2 stand for the effective 
theory 



cff 



1 

2m 



2mc 2 



d ip ■ 



1 

2m 



(3.22) 



where we recall that n is the particle number density of the fluid at rest. The theory 
describes a nonrelativistic sound wave, with the dimensionless phase field tp represent- 
ing the Goldstone mode of the spontaneously broken global U(l) symmetry. It has the 
gapless dispersion relation E 2 {k) = c 2 k 2 . The effective theory gives a complete de- 
scription of the superfluid valid at low energies and small momenta. The same effective 
theory appears in the context of (neutral) superconductors [|l3[| (see next chapter) and 
also in that of classical hydrodynamics 

The chemical potential /i is represented in the effective theory ( 3.22 ) by [ |l5t 



fi(x) = -d tp(x), 



so that 



d£ e s d£ cS 



Of l 



n(x), 



(3.23) 



(3.24) 
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as required. It also follows from this equation that the particle number density n(x) is 
canonical conjugate to — <f>(x). 



The most remarkable aspect of the effective theory ( 3.22 ) is that it is nonlinear. The 
nonlinearity is necessary to provide a Galilei-invariant description of a gapless mode, 
as required in a nonrelativistic context. Under a Galilei boosts, 

t^t' = t, x^x'=x-ui; d -> d' = d + u- V, V -> V' = V, (3.25) 

with u a constant velocity, the Goldstone field <p{x) transforms as 

—<p(x) -> —<p'(x') = —<p(x) - u ■ x + ±u 2 i. (3.26) 

m m m z 

As a result, the superfluid velocity v s = Wip/rn and the chemical potential (per unit 
mass) /i/m = ~doip/m transform under a Galilei boost in the correct way, 

v s (x) — > Vg(x') = v s (a;) — u, //(x)/m — > f/(x')/m = (i{x)/m — u ■ v s (x) + ^u 2 . 

(3.27) 

It is readily checked that the field U(x) defined in ( 3.10[ ) and therefore the effective 



theory (3.22) is invariant under Galilei boosts. 



Since the Goldstone field in (3. 1 ) is always accompanied by a derivative, we see that 
that the nonlinear terms carry additional factors of |k| /mc, with |k| the wave number. 
They can therefore be ignored provided the wave number is smaller than the inverse 
coherence length = mc, 

|k| < 1/e (3.28) 
For example, in the case of 4 He the coherence length, or Compton wavelength, is 



about 10 nm. In this system, the bound (3.28 ), below which the nonlinear terms can be 



neglected, coincide with the region where the spectrum is linear and the description in 
terms of solely a sound mode is applicable. 

The alert reader might be worrying about an apparent mismatch in the number of 
degrees of freedom in the normal and the superfluid phase. Whereas the normal phase 
is described by a complex 0-field, the superfluid phase is described by a real scalar 
field ip. The resolution of this paradox lies in the spectrum of the modes Jl^]. In 
the normal phase, the spectrum E(\s.) = k 2 /2m is linear in E, so that only positive 
energies appear in the Fourier decomposition, and one needs — as is well known from 
standard quantum mechanics — a complex field to describe a single particle. In the 
superfluid phase, where the spectrum, E 2 (lt) = c 2 k 2 , is quadratic in E, the counting 
goes differently. The Fourier decomposition now contains positive as well as negative 
energies and a single real field suffice to describe this mode. In other words, although 
the number of fields is different, the number of degrees of freedom is the same in both 
phases. 



The particle number density and current that follows from (3.22) read 



n ( g ^ 1 
mc 2 \ 2m 1 

j(x) = n(x)v s (x). (3.30) 



"(*) = ^-Z^\^(x) + —[V l p(x)} 2 \ (3.29) 
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Physically, (3.29) reflects Bernoulli's principle which states that in regions of rapid 



flow, the density and therefore the pressure is low. 



The diagrams of Fig. 3.2 can be evaluated in a loop expansion to obtain explicit ex- 



pressions for the particle number density n and the sound velocity c to any given order 
fll2[]. In doing so, one encounters — apart from ultraviolet divergences which will be 
dealt with shortly — also infrared divergences because the Bogoliubov spectrum is gap- 
less. When however all one-loop contributions are added together, these divergences 
are seen to cancel [jO]. One finds for d — 2 to the one-loop order 



2 A' 



An 



(3.31) 



where u and A are the renormalized parameters. Following Ref. Jl7[], we adopted a 
dimensional regularization scheme, in which after the integrals over the loop energies 
have been carried out, the remaining integrals over the loop momenta are analytically 
continued to arbitrary space dimensions d. As renormalization prescription we em- 
ployed the modified minimal subtraction scheme. This leads to the following relation 



between the bare (Aq) and renormalized coupling constant [see Eq. (3.53) below] 



1 



1 (i 



rn 



(3.32) 



where e — 2 — d, and n is an arbitrary renormalization group scale parameter introduced 
to give the renormalized coupling constant A the same engineering dimension as in 
d = 2. The chemical potential is not renormalized to this order. 

Incidentally, from the vantage point of renormalization, the mass m is an irrelevant 
parameter in nonrelativistic theories and can be scaled away (see, e.g., Ref. [|l7||). 

The form of the effective theory (3.22) can also be derived from general symmetry 
arguments [|l8|], More specifically, it follows from making the presence of a gapless 
Goldstone mode compatible with Galilei invariance which demands that the mass cur- 
rent and the momentum density are equal. The latter observation leads to the conclu- 
sion that the U(l) Goldstone field tp can only appear in the combination (3.10). To 
obtain the required linear spectrum for the Goldstone mode it is necessary then to have 
the form (3.22). Given the form of the effective theory, the particle number density 



and sound velocity can then more easily be obtained directly from the thermodynamic 
potential Q, via 

i an i i md 2 n 

V o/j, cr V n op? 

where V is the volume of the system. In this approach, one only has to calculate the 
thermodynamic potential which at zero temperature and in the Bogoliubov approxima- 
tion in which we are working is given by the sum V of the classical potential Vo and 
the effective potential V c g corresponding to the theory (3.1): 



(Vo + Veff) 



(3.34) 



where Vo is given by (3.4 ) with <j) replaced by <f>. The effective potential for the uniform 
system is obtained as follows. In the Bogoliubov approximation of ignoring higher than 
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second order in the fields, the integration over $ is Gaussian. Carrying out this integral, 
we obtain for the zero-temperature partition function 



' J> J D0*D0exp [ij^ {2) 
^DeiT^Mo), 



where Mq stands for the matrix introduced in (p\9f). Setting 



Z = cxp 



i - / V + S, 



oil' 



(3.35) 



(3.36) 



we conclude from (3.35) that the effective action in the Bogoliubov approximation is 

y 

Seff = iiTrm[M (p,ar)], (3.37) 



given to the one-loop order by 



where we again used the identity Det(A) = exp[Tr ln(A)]. The trace Tr appearing 
here stands besides for the trace over discrete indices now also for the integral J over 
spacetime as well as the one J k over energy and momentum. The latter integral reflects 
the fact that the effective action calculated here is a one-loop result with k^ the loop 
energy and momentum. To disentangle the integrals one has to carry out similar steps 
as the ones outlined in Sec. 12 and repeatedly apply the identity 



f{x)pfj,g(x) = - id ll )f(x)g(x), 



(3.38) 



where f(x) and g[x) are arbitrary functions of spacetime and the derivative <9 M = 
(do, —V) acts only on the next object to the right. The method outlined there can 
easily be transcribed to the present case where the time dimension is included. 

If the field U(x) in Mo is set to zero, things simplify because Mo now depends on 
Pfj, only. The effective action then becomes S c s = — f V c e with 



I 

--tr 
2 



ln[M (fc)] 



(3.39) 



the effective potential. The easiest way to evaluate the integral over the loop variable k 
is to first differentiate the expression with respect to the chemical potential /ig 



— tr / ln[M (A)] = -2 



e(k) 



£ 2 (k) 



(3.40) 



with E(k) the Bogoliubov spectrum (3.13). The integral over fco can be carried out 
with the help of a contour integration, yielding 



5(k) 



k k 2 -E*(k) 



17] 



1 f £ ( k ) 

2 Aw 



(3.41) 



This in turn is easily integrated with respect to /i . Putting the pieces together, we 
obtain 

1 



V 



(3.42) 
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The integral over the loop momentum in arbitrary space dimension d yields 

T(l - d/2)r(d/2 + l/2) 



V 



9 

M5 



T d/2 d/2+1 

L d m 1 fi ' 



(3.43) 



4A """" r ' u ' " d 27T d / 2 + 1 /2r(d/2 + 2) 
where we employed the integral representation of the Gamma function 

i = j_ r *> e -~ 

a z TO) y r 

together with dimensional regularization to suppress irrelevant ultraviolet divergences. 

For comparison, let us also evaluate the integral in (3.42) over the loop momentum 
in three dimensions by introducing a momentum cutoff A 



(3.44) 



V, 



i 

47T 2 



1A 5 1 , 3 2 a 32 

— h -MoA - m/i A + —to 

10 to 3 15 



?'^T)+o{^. (3.45) 



From ( |3.43| ), we obtain by setting d = 3 only the finite part, so that all terms diverging 
with a s trict ly positive power of the momentum cutoff are suppressed. As we remarked 
in Sec. 2.2, these contributions, which come from the ultraviolet region, cannot be 
physically very relevant because the simple model (3.1) breaks down here. On account 
of the uncertainty principle, stating that large momenta correspond to small distances, 
these terms are always local and can be absorbed by redefining the parameters appear- 
ing irrthe Lagrangian [JL9]]. Since = 2Ao|0| 2 , we see that the first diverging term 
in (3.45) is an irrelevant constant, while the two remaining diverging terms can be 
absorbed by introducing the renormalized parameters 



1 „ 

A = A 2 toAqA. 



(3.46) 
(3.47) 



Because the diverging terms are — at least to this order — of a form already present in the 
original Lagrangian, the theory is called "renormalizable". The renormalized parame- 
ters are the physical ones that are to be identified with those measured in experiment. 
In this way, we see that the contributions to the loop integral stemming from the ultra- 
violet region are of no importance. What remains is the finite part 



V eff = 



157T 



™3/2 5/2 



(3.48) 



which, as we have seen, is obtained directly without renormalization when using di- 
mensional regularization. In this scheme, divergences proportional to powers of the 
cutoff never show up. Only logarithmic divergences appear as 1/e poles, where e is the 
deviation from the upper critical dimension (d = 2 in the present case). These loga- 
rithmic divergences \n(A/E), with E an energy scale, are relevant also in the infrared 
because for fixed cutoff \n(A/E) — > — oo when E is taken to zero. 
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In so-called "nonrenormalizable" theories, the ultraviolet-diverging terms are still 
local but not of a form present in the original Lagrangian. Whereas in former days 
such theories were rejected because there supposed lack of predictive power, the mod- 
ern view is that there are no fundamental theories and that there is no basic difference 
between renormalizable and nonrenormalizable theories [ p0| ] . Even a renormalizable 
theory like (3.1) should be extended to include all higher-order terms such as a |0| 6 - 



term which are allowed by symmetry. These additional terms render the theory "non- 
renormalizable". This does not however change the predictive power of the theory. 
The point is that when describing the physics at an energy scale E far below the cutoff, 
the higher-order terms are suppressed by powers of E/A, as follows from dimensional 
analysis. Therefore, far below the cutoff, the nonrenormalizable terms are negligible. 

That d = 2 is the up per critical dimension of the problem at hand can be seen by 
noting that L d in ( 3.43 ) diverges when d tends to 2. Special care has to be taken for 
this case. For d ^ 2, we obtain with the help of ( 3.33 ) [ pT] ] 

Mo 



2A 



l + (d + 2)L d m d / 2 \ of i?/ 2 - 1 



and 



Mo 
rn 



1 



{d-2){d/2 + l)L d m d l 2 \ m d J 2 - 



(3.49) 



(3.50) 



where to arrive at the last equation an expansion in the coupling constant Ao is made. 
Up to this point, we have considered the chemical potential to be the independent pa- 
rameter, thereby assuming the presence of a reservoir that can freely exchange particles 
with the system under study. The system can thus have any number of particles, only 
the average number is fixed by external conditions. From the experimental point of 
view it is, however, often more realistic to consider the particle number fixed. If this is 
the case, the particle number density n should be considered as independent variable 
and the chemical potential should be expressed in terms of it. This can be achieved by 



inverting relation (3.49): 



Ho = 2A n [l - 2(d - 2){d/2 + l)L d m d/2 \ Q (2\ n) d/2 
The sound velocity expressed in terms of the particle number density reads 

2Xofl [ 1 - d(d/2 + l)L d m d ' 2 X Q (2\ a n) d l 2 - 



These formulas reproduce the known results in d = 3 [H] and d = 1 



(3.51) 



(3.52) 



To investigate the case d = 2, we expand the potential ( |3.43| ) around d = 2: 

4A 47T6 K e 

with e = 2 — d. This expression is seen to diverge in the limit d — ► 2. The theory 



(3.53) 



can be rendered finite by introducing a renormalized coupling constant via ( 3.32 ). We 
also see that the chemical potential is not renormalized to this order. The beta function 
/3(A) follows as 



-eA H A 

7T 



(3.54) 
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In the upper critical dimension, this yields only one fixed point, viz. the infrared-stable 
(IR) fixed point A* = 0. Below d = 2, this point is shifted to A = eir/m. It is now 



(3.31) 



easily checked that Eqs. (3.49) and (3.52) also reproduce the two-dimensional results 



In the one-loop approximation there is no field renormalization; this is the reason 



why in (3. 1 ) we gave only the bare parameters /io and Aq an index 0, and not > 



We proceed by calculating the fraction of particles residing in the condensate, 
deriving the Bogoliubov spectrum (3.13), we set \cf>\ 2 



number density of particles contained in the condensate, 

no = H 2 , 



In 



Ho/2Xq thereby fixing the 



(3.55) 



in terms of the chemical potential. For our present consideration we have to keep <fi as 
independent variable. The spectrum of the elementary excitation expressed in terms of 

d> is 



E(k) = y/ [e(k) - Mo + 4AqH 2 ] - 4A§|0|4 . (3.56) 

It reduces to the Bogoliubov spectrum when the mean-field value (|3.6|) for </> is inserted. 
Equation (3.39) for the effective potential is still valid, and so is (3.34). We thus obtain 
for the particle number density 



1 d 

2 dfi 



E(k) 



(3.57) 



0l 2 =A t o/2A o 



where the mean-field value for <f> is to be substituted after the differentiation with re- 
spect to the chemical potential has been carried out. We find 



^2 _ 2^/2-2^^/2^/2^ 

a — 1 



(3.58) 



or for the so-called depletion of the condensate [E5p 



n 
no 



1 



_2<i/2 



1 



■L d m d l 2 \ d / 2 n d / 2 -\ 



(3.59) 



where in the last term we replaced the bare coupling constant with the (one-loop) renor- 
malized o ne. Th is is consistent to this order since this term is already a one-loop result. 
Equation ( 3.59 ) shows that even at zero temperature not all the particles reside in the 
condensate. Due to the interparticle repulsion, particles are removed from the zero- 
momentum ground state and put in states of finite momentum. It has been estimated 
that in bulk superfluid 4 He — a strongly interacting system — only about 8% of the par- 
ticles condense in the zero-momentum state [ 26 1 . For d = 2, the right-hand side of Eq. 
( 3.59 ) reduces to 

^ - 1 « (3.60) 

no Zir 

which is seen to be independent of the particle number density. 

Despite the fact that not all the particles reside in the condensate, they all participate 
in the superfluid motion at zero temperature [E7|. Apparently, the condensate drags 
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the normal fluid along with it. To show this, let us assume that the entire system 
moves with a velocity u relative to the laboratory system. As is known from standard 
hydrodynamics the time derivate in the frame following the motion of the fluid is 8q + 



u • V [see Eq. (3.25)]. If we insert this in the Lagrangian (3.1) of the interacting Bose 
gas, it becomes 

C = 0*[id o - e(-tV) + Mo - u • (-»V)]0 - A o |0| 4 , (3.61) 

where the extra term features the total momentum f </>*(— zV)</> of the system. The 
velocity u multiplying this is on the same footing as the chemical potential Mo mul- 
tiplying the particle number J x \tp\ 2 . Whereas n is associated with particle number 
conservation, u is related to the conservation of momentum. 

In the two-fluid picture, the condensate can move with a different velocity v s as the 



rest of the system. To bring this out we introduce new fields, cf. (3.5) 

4>{x) (j)'(x) = e lmv "- x (t)(x) (3.62) 

in terms of which the Lagrangian becomes fl28| ] 

C = 0* [id - e(-*V) +Mo - \rnw s • (v s - 2u) - (u - v.) ■ (-<V)] tj>- A o |0| 4 , (3.63) 

where we dropped the primes on <fi again. Both velocities appear in this expression. 
Apart from the change u — > u — v s in the second last term, the field transformation 
resulted in a change of the chemical potential 

Mo — ► Mcff := Mo - ^rnvg • (v s - 2u) (3.64) 

where Meff may be considered as an effective chemical potential. 

The equations for the Bogoliubov spectrum and the thermodynamic potential are 
readily written down for the present case with these two changes are kept in mind. In 



particular, the effective potential is given by (3.42) with the replacement Eq. (3.64). 
The momentum density, or equivalently, the mass current g of the system is obtained 
in this approximation by differentiating the effective potential with respect to — u. We 
find, using the equation 

— — = mv s (3.65) 
on 

that it is given by 

g = A=v s , (3.66) 

with p s = mn the superfluid mass density. This equation, comprising the total particle 
number density n, shows that at zero temperature indeed all the particles are involved 
in the superflow, despite the fact that only a fraction of them resides in the condensate 
|^7|]. The superfluid mass density p s , obtained by evaluating the response of the system 
to an externally imposed velocity field u, should not be confused with the number 



density no of particles contained in the condensate introduced in Eq. (3.55) 



Let us close this section by pointing out a quick trail to arrive at the effective theory 



(3.22) starting from the microscopic model (3.1). To this end we set 



0(a) =e iv ^[Vn+4>(x)}, (3.67) 
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and expand the Lagrangian (3.1) up to quadratic terms in <j). This leads to 

£( 2 ) = -V -hU- VEU(4> + <j>*) - X Q n{4> + 4>*f, (3.68) 

where we used the mean-field equation /j,q = 2Aon. We next integrate out the tilde 
fields — which is tantamount to substituting the field equation for these fields back into 
the Lagrangian — to obtain 

C cS = -fiU{x) + -U(x) — U(x), (3.69) 

4 An 

apart from the irrelevant constant term Vo- This form of the effective theory is equiv- 
alent to the one found before in (3.1 ). We have cast the last term in a form that can 
be easily generalized to systems with long-ranged interactions. A case of particular 
interest to us is the Coulomb potential 

V(x) = ^r, (3.70) 

whose Fourier transform in d space dimensions reads 

V(k) = 2 d -V rf - 1 >/ 2 r [i(d- 1)] ppg^. (3.71) 



The simple contact interaction Li = — Xq J x |</>(a;)| 4 in (3.1) gets now replaced by 

Li = ~ f m^)\ 2 V( X -y)m,y)\ 2 . (3.72) 

The rationale for using the three-dimensional Coulomb potential even when consider- 
ing charges confined to move in a lower dimensional space i s tha t the electromagnetic 
interaction remains three-dimensional. The effective theory ( 3.69 ) now becomes in the 
Fourier representation 

£ cS = -nU(k) + ±U(k ,-k)^-U(k , -k) (3.73) 

and leads to the dispersion relation 

E 2 (k) = 2V d - 1 )/ 2 r [I(d - 1)1 ^|k| 3 ~ d . (3.74) 

For d = 3, this yields the famous plasma mode with an energy gap given by the plasma 
frequency Wp = Airne^/m. 

To appreciate under which circumstances the Coulomb interaction becomes impor- 
tant, we note that for electronic systems l/|x| ~ fcp for dimensional reasons and the 
fermion number density n ~ hp, where fcp is the Fermi momentum. The ratio of the 
Coulomb interaction energy ec to the Fermi energy ep = fcj?/2m is therefore pro- 
portional to nT x l d . This means that the lower the electron number density, the more 
important the Coulomb interaction becomes. 
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3.3 Quenched Impurities 

In most of the quantum systems we will be considering, impurities plays an impor- 
tant role. The main effect of impurities is typically to localize states. Localization 
counteracts the tendency of the system to become superfluid. We shall therefore now 
include impurities in the interacting Bose gas to see whether this leads to localization 
and whether the system still has a superfluid phase. It is expected that on increasing the 
strength of the disorder for a given repulsive interparticle interaction, the superfluid un- 
dergoes a zero-temperature phase transition to an insulating phase of localized states. 
The location and nature of this transition will be the subject of Ch. [| 

We shall assume that the impurities are fixed and that their distribution is not af- 
fected by the host system. This type of impurities is called quenched impurities and is 
to be distinguished from so-called annealed impurities which change with and depend 



on the host system. To account for impurities, we add to the theory (3.1) the term 

£ A =^(x)|0(x)| 2 , (3.75) 
with ^>(x) a random field whose distribution is assumed to be Gaussian [E9| 



P(V>) = cxp 



(3.76) 



and characterized by the disorder strength Ao. The engineering dimension of the ran- 
dom field is the same as that of the chemical potential which is one, [tp] = 1, while that 



of the parameter Ao is [Ao] = 2 — d so that the exponent in (3.76) is dimensionless. 
Since ^(x) depends only on the d spatial dimensions, the impurities it describes should 
be considered as grains randomly distributed in space. The quantity 



Zty] = / D0*I 



D(j) exp [i / £ , (3.77) 



where now C stands for the Lagrangian (3.1 ) with the term (3.75) added, is the zero 



temperature partition function for a given impurity configuration ip. In the case of 
quenched impurities, the average of an observable 0(0*, 4>) is obtained as follows 

(0(0*, 4>)) = [ DV>PWOW*,0)>^ (3.78) 



where (0(<fr* , (f)) ^ indicates the grand-canonical average for a given impurity config- 
uration. In other words, first the ensemble average is taken, and only after that the 
averaging over the random field is carried out. 

In terms of the shifted field, the added term reads 

£ A = V>(x)(M 2 + |0| 2 + H>* + P4>)- (3.79) 

The first two terms lead to an irrelevant change in the chemical potential, so that we 
only have to consider the last two terms, which we can cast in the form 

C A = i/>(x)&§, $= ( t V (3.80) 
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The integral over $ is Gaussian in the Bogoliubov approximation and is easily per- 
formed to yield an additional term to the effective action 



Sa = 



(3.81) 



:r,y 



where the propagator Go is the inverse of the matrix M introduced in (3.9) with the 
field U (x) set to zero. Let us first Fourier transform the fields, 



e -ik-(x-y) Q Q ( k j 

e 4k - x ^(k). 



Go(x -y) = 
V-(x) = 

The contribution to the effective action then appears in the form 

SA = -^|<Mk)| 2 i>tG(o,k)i>. 



(3.82) 
(3.83) 

(3.84) 



Since the random field is Gaussian distributed [see (3.76)], the average over this field 
representing quenched impurities yields, 



(hA(k)| 2 ) = iT/A„. 



(3.85) 



The remaining integral over the loop momentum in (3.84) is readily carried out to yield 



1 



,\d/2 



(£ A > = -T(l - d/2) j |0| 2 (6AoH 

This contribution is seen to diverge in the limit d — > 2: 

1 m/io A 



Mo 



)<*/2-l Ao . 



47T AoK e e 



(3.86) 



(3.87) 



where we substituted the mean-field value /io = 2Ao|</>| 2 . Recall that k is an arbitrary 
scale parameter introduced for dimensional reasons; the engineering dimension of the 
right-hand side in (3.87) has the correct value 3 — e this way. The result (3.87) is a first 
indication of the importance of impurities in d = 2, showing that in order to render the 
random theory finite a modified renormalized coupling constant A has to be introduced 
via, cf. (|3.32|), 

A ~ 



1 
Ao 



1 

K' 




(3.88) 



which depends on the disorder strength. The renormalized parameter A is defined in 
the same way as A. 

In the previous section we saw that due to the interparticle repulsion, not all the 
particles reside in the condensate. We expect that the random field c auses an additional 
depletion of the condensate. To obtain this, we differentiate ( 3.86 ) with respect to the 
chemical potential. This gives [ 30 1 



dfi 



2d/2 " T{2 ~ d,2 ± m d / 2 \ d l 2 -*n d ' 2 - 1 & 1 (3.89) 



T d/2 
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where no denotes the density of particles residing in the condensate. We have here 
again replaced the bare parameters with the (one-loop) renormalized ones. This is 



consistent to this order since (3.89 ) is already a one-loop result. 

The divergence in the limit A — > for d < 4 signals the collapse of the system when 
the interparticle repulsion is removed. Note that in d = 2, the depletion is independent 



of the condensate density no [31 1: 



1 771 
167T A 

The total particle number density n is given by 

/ mA\ 1 rn . 
n — n (l + — — j + — — —A. (3.91) 

We next calculate the mass current g to determine the superfluid mass density, 
i.e., the mass density flowing with the superfluid velocity v s . As we have seen in the 
preceding section, in the absence of impurities and at zero temperature all the particles 
participate in the superfiow and move on the average with the velocity v s . We expect 
this to no longer hold in the presence of impurities. To determine the change in the 



superfluid mass density due to impurities, we replace /jq with u e g as defined in (3.64) 



and ido with ido — (u — v s ) • (— iV) in the contribution ( 3.84 ) to the effective action, 
and differentiate it with respect to — u — the externally imposed velocity. We find to 
linear order in the difference u v s : 

g = p s v s + p n u, (3.92) 



with the superfluid and normal mass density [ 30 1 



p s = m - ^ n ^j ' P a = ^mn A . (3.93) 

We see that the normal density is a factor A/d larger than the mass density mn& 
knocked out of the condensate by the impurities. (For d = 3 this gives the factor | 
first found in Ref. p2|].) Apparently, part of the zero-momentum states belongs for 
d < 4 not to the condensate, but to the normal fluid. Being trapped by the impuri- 
ties, this fraction of the zero-momentum states are localized. This shows that the phe- 
nomenon of localization can be accounted for in the Bogoliubov theory of superfluidity 
by including a random field. 



3.4 Vortices 

We shall now include vortices in the system. A vortex in two space dimensions may 
be pictured as a point-like object at scales large compared to their core seize. It is 
characterized by the winding number w of the map 

(p(x) : Si -» S 1 (3.94) 
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of a circle around the vortex into the internal circle S 1 parameterized by the Gold- 



stone field ip. In the microscopic theory (3.1 ), the asymptotic solution of a static vortex 



with winding number w located at the origin is well known [ 33 1 



2A V 4x2 / V x 

where 9 is the azimuthal angle and £o = 1/ \/™4k> = 1 / mc o is the coherence length. 
The density profile n(x) in the presence of this vortex follows from taking |</>(x) | 2 . 

To incorporate vortices in the effective theory we employ the powerful principle 
of defect gauge symmetry developed by Kleinert [Q, [36|. In this approac h, one 
introduces a so-called vortex gauge field tp p = ((p 9 , <p p ) in the effective theory ( 3.22j ) 
via minimally coupling to the Goldstone field: 

dp<p -> dn<P + <p*, (3-96) 

with 9^ = (do, —V). If there are N vortices with winding number w a (a = 1, • • • , N) 
centered at X 1 (*),•••, ~K N (t), the plastic field satisfies the relation 

V x ip p (x) = -2vr^w Q( 5[x~X a (i)], (3.97) 

a 

so that we obtain for the superfluid velocity field 

Vxv s = ^ 7ll %-r(i)], (3.98) 

a 

as required. Here, 7 Q = (2tt /m)w a is the circulation of the ath vortex which is 
quantized in units of 27r/m. A summation over the indices labeling the vortices will 
always be made explicit. The combination d^ip + ip p is invariant under the local gauge 
transformation 

ip(x) -> ip(x) + a(x); ifl^ipl~ d^a{x), (3.99) 
with (p p playing the role of a gauge field. 



In the gauge <p^ = 0, Eq. (3.97) can be solved to yield 



tf(x) = 2Tre ij ^2w a S j [x,L a (t)] (3.100) 

a 

where e»j is the antisymmetric Levi-Civita symbol in two dimensions, with e\2 = 1, 
and 6[x, L a (t)] is a delta function on the line L a (t) starting at the center X a (t) of the 
ath vortex and running to spatial infinity along an arbitrary path: 

8i[x,L a {t)} = f di/i«(x-y). (3.101) 

JL a (t) 

Let us for the moment concentrate on static vortices. The field equation obtained 



from the effective theory (3.22) with \Jtp replaced by the covariant derivative VtyS — cp 



p 



and dap set to zero simply reads 

V • v s = 0, or V • (Vip - <p p ) = 0, (3.102) 
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when the fourth-order term is neglected. It can be easily solved to yield 



p(x) = - / G(x-y)V- V p (y), 



where G(x) is the Green function of the Laplace operator 

2~n 



G(x) = 



For the velocity field we obtain in this way the well-known expression [B7p 



1 N x- 
Wi (x) = - £ij .^ 7a -2 

Q = l ' 



■x? 



-X Q | 2 ' 



(3.103) 



(3.104) 



(3.105) 



which is valid for x sufficiently far away from the vortex cores. Let us now specialize 
to the case of a s ingle static vortex at the origin. On substituting the corresponding 
solution in ( 3.29 ), we find for the density profile in the presence of a static vortex 
asymptotically 



(3.106) 



This is the same formula as the one obtained from the solution ( 3.95 ) of the microscopic 



theory. This exemplifies that with the aid of the defect gauge symmetry principle, 
vortices are correctly accounted for in the effective theory. 

Let us proceed to investigate the dynamics of vortices in this formalism and de- 
rive the action which governs it. We consider only the first part of the effective theory 



(3.22). In ignoring the higher-order terms, we approximate the superfiuid by an in- 
compressible fluid for which the particle number density is constant, n(x) = n, see 
Eq. (3.29). We again work in the gauge ipg = and replace Wip by the covariant 
derivative Vip~ ip p , with the plastic field given by ( 3.97[ ). The solution of the resulting 
field equation for ip is again of the form ( 3.103 ), but now it is time-dependent because 
the plastic field is. Substituting this in the action S e g — f £ e ff> we find a ft er some 
straightforward calculus 



S e B = 



(3.107) 



The constant a has the dimension of a length and is included in the argument of the 
logarithm for di mension al reasons. Physically, it represents the core size of a vortex. 
The first term in (3.107) leads to a twisted canonical structure which is reminiscent of 



that found in the so-called Landau problem of a charged particle confined to move in a 
plane perpendicular to an applied magnetic field H. 

To display the canon ical structure, let us rewrite the first term of the Lagrangian 
corresponding to (3.107) as 



Li =mnV 7a Xfi 2 a , 



(3.108) 
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where we ignored a total derivative. It follows that the canonical conjugate to the 
second component X£ of the center coordinate X Q is essentially its first component 
PI 

9ii 



dX? 



mnj a X?. (3.109) 



It implies that phase space coincides with real space and gives rise to the commutation 
relation 

[x 2 Q ,xf] = ^-e 2 s at3 , (3.iio) 

where 

e=l/V2nh (3.111) 

is a characteristic length whose definition is such that 2tt£ 2 is the average area occupied 
by a particle of the superfluid film. The commutation relation leads to an uncertainty 
in the location of the vortex centers given by 

AX? AX? > ^— (3.112) 

2\w a \ 

which is inverse proportional to the particle number density. 

From elementary quantum mechanics we know that to each unit cell (of area h) in 
phase space there corresponds one quantum state. That is, the number of states in an 
area S of phase space is given by 

#statesin S = \ [ dpdq, (3.113) 
hj s 

where p and q are a pair of canonically conjugate variables. For the case at hand, this 
implies that the available number of states in an area S a of real space is 

# states in S a = \w a \ nS a , (3.114) 

or, equivalently, that the number of states per unit area available to the ath vortex is 
\w a \n. 

This phenomenon that phase space coincides with real space is known to also arise 
in the Landau problem. There, it leads to the well-known degeneracy \e a \ H/h of each 
Landau level, where e a = v a eo is the electric charge of the particle, with e (> 0) 
the unit of charge. In terms of the magnetic flux quantum <f> = h/e , the Landau 
degeneracy can be rewritten as \v a \H/ $o = I^q with the flux number density. 
In other words, whereas the degeneracy in the case of vortices in a superfluid film 
is given by the particle number density, here it is given by the flux number density. 
Using this analogy, we see that the characteristic length ( |3.111[ ) translates into £h = 
1 / y / 2~Kn^ which is precisely the magnetic length of the Landau problem. 



The first term in the action (3.107) is also responsible for the so-called geometrical 
phase [ p9| ] acquired by the wavefunction of a vortex when it traverses a closed path. 
Let us first discuss the case of a charged particle moving adiabatically around a close 
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path T Q . Its wavefunction picks up an extra Aharonov-Bohm phase factor given by the 
Wilson loop: 



W(T a ) = exp[i 7 (r Q )] = exp ( ^ 4 dx • A 



exp 



2mv, 



HS(T a ) 



(3.115) 



where A is the vector potential describing the external magnetic field and HS(T a ) is 
the magnetic flux through the area S(T a ) spanned by the loop T a . The geometrical 
phase 7 (r Q ) in (3.115) is seen to be (2Trv a times) the number of flux quanta enclosed 
by the path T a . 

On account of the above analogy, it follows that the geometrical phase picked up 
by the wavefunction of a vortex when it is moved adiabatically around a closed path in 
the superfluid film is (2irw a times) the number of superfluid particles enclosed by the 
path [@. 



The second term in the action (3.107) represents the long-ranged interaction be- 
tween two vortices mediated by the exchange of Goldstone quanta. The action yields 
the well-known equations of motion for point vortices in an incompressible two-dimensional 
superfluid @fu|]: 



*f (*) = r ^ 



V- *f(*) 
7a 



|X0(t)-X a (t)| s 



(3.116) 



Note that X?(t) = [X /3 (t)] , where v(x) is the superfluid velocity ( |3.105[ ) with the 
time-dependence of the centers of the vortices included. This nicely illustrates a result 
due to Helmholtz for ideal fluids, stating that a vortex moves with the fluid, i.e., at the 
local velocity produced by the other vortices in the system. Experimental support for 
this conclusion has been reported in Ref . [B2F . 



3.5 Kosterlitz-Thouless Phase Transition 

Although we are interested mainly in quantum phase transitions in these Lectures, there 
is one classical phase transition special to two dimensions which turns out to be relevant 
for our discussion later on — the so-called Kosterlitz-Thouless phase transition. It is 
well known that a superfluid film undergoes such a phase transition at a temperature 
well below the bulk transition temperature. The superfluid low-temperature state is 
characterized by tightly bound vortex-antivortex pairs which at the Kosterlitz-Thouless 
temperature unbind and thereby disorder the superfluid state. The disordered state, 
at temperatures still below the bulk transition temperature, consists of a plasma of 
unbound vortices. 

Since the phase transition is an equilibrium transition, we can ignore any time de- 
pendence. The important fluctuations here, at temperatures below the bulk transition 
temperature, are phase fluctuations so that we can consider the London limit, where the 
phase of the (/)(x)-field is allowed to vary in spacetime while the modulus is kept fixed, 
and take as Hamiltonian 

H=± Ps vl (3.117) 
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where p s is the superfluid mass density which we assume to be constant and v s is the 
superfluid velocity 



(3.118) 



with the vortex gauge field (p p included to account for possible vortices in the system. 
We shall restrict ourselves to vortices of unit winding number, so that w a = ±1 for a 
vortex and antivortex, respectively. 

The canonical partition function describing the equilibrium configuration of N + 
vortices and i\T_ antivortices in a superfluid film is given by 



1 



iV + !AU 



n 



cxp 1-/3 H 



(3.119) 



with TL the Hamiltonian ( 3.117| ) and N = N + + N_ the total number of vortices 
and antivortices. The factors N+\ and N-\ arise because the vortices and antivortices 
are indistinguishable, and J\ a L. a denotes the integration over the positions of the 
vortices. The functional integral over ip is Gaussian and therefore easily carried out, 
with the result 



1 



N+\N_\ 



n 



exp 



7r %£«, a «, /J ]n(|X«-X'»|/a) 



a,/3 



(3.120) 



Apart from an irrelevant normalization factor, Eq. ( |3.120| ) is the canonical partition 
function of a two-dimensional Coulomb gas with charges q a = qw a = ±q, where 



(3.121) 



Let us rewrite the sum in the exponent appearing in (3.120) as 



1 a qp [In (|X« - X^|/ a ) - ln(0)] + ln(0) ^ q a 



(3.122) 



where we isolated the self-interaction in the last term at the right-hand side. Since 
ln(0) = — oo, the charges must add up to zero so as to obtain a nonzero partition 
function. From now on we will therefore assume overall charge neutrality, ^ Q q a = 0, 
so that N + = N- = N/2, where N must be an even integer. To regularize the 
remaining divergence, we replace ln(0) with an undetermined, negative constant — c. 
The exponent of ( 3.120| ) thus becomes 

| QaQp In (|X Q - X"|/a) = f £ q a q P In (|X« - X*|/a) - ?e c N, (3.123) 



a, /3 



a^f3 



where e c = cq 2 /2 physically represents the core energy, i.e., the energy required to 
create a single vortex. In deriving this we used the identity J2 a ^/3 QaQp = ~ J2 a Qa = 
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-Nq 2 which follows from charge neutrality. Having dealt with the self-interaction, 



we limit the integrations J\ a J Xa in (3.12C) over the location of the vortices to those 
regions where they are more than a distance a apart, |X a — X^| > a. The grand- 
canonical partition function of the system can now be cast in the form 



z = 



oo 



n 



exp 







(3.124) 



where z — exp(— j3e c ) is the fugacity. The system is known to undergo a phase transi 
tion at the Kosterlitz-Thouless temperature 

1 . 



;ystem is 



KT 



7T p s 

2 m 2 ' 



(3.125) 



triggered by the unbinding of vortex-antivortex pairs. It follows from this equation 
that the two-dimensional superfluid mass density p s {T), which varies from sample to 
sample, terminates on a line with universal slope as the temperature approaches the 
Kosterlitz-Thouless temperature from below Mq]. 



3.6 Dual Theory 



Let us proceed to represent the partition function ( 3 . 1 24 ) by a field theory — a so-called 



dual theory. The idea behind such a dual transformation is to obtain a formulation in 
which the vortices are not described as singular objects as is the case in the original 
formulation, but by ordinary fields. To derive it we note that ln(|x|) is the inverse of 
the Laplace operator V 2 , 

^V 2 ln(|x|) = <5(x). (3.126) 



This allows us to represent the exponential function in ( 3.124 ) as a functional integral 
over an auxiliary field </>: 



exp 



^ fo(Z/3 ln(|X Q -X^A 



D</>exp < — 



1 

4tt/3 



(V0) 2 + i Pq <i> 



(3.127) 

where p 9 (x) = J2 a — X Q ) is the charge density. In this way, the partition 

function becomes 



N 

AJ=0 LV ' ' 1 a = l 



Dtp exp < — 



(3.128) 



In a me an-field treatment, the functional integral over the auxiliary field introduced in 
(3.127) is approximated by the saddle point determined by the field equation 



-2itp q . 



(3.129) 
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When we introduce the scalar variable $ := iTip, this equation becomes formally 
Gauss' law, with <i> the electrostatic scalar potential. The auxiliary field introduced in 



(3.127) may therefore be thought of as representing the scalar potential of the equiva- 



lent two-dimensional Coulomb gas [|34 

On account of charge neutrality, we have the identity 



N 



Nl 



{(N/2) 



M2 



N 

n 



(3.130) 



where we recall that TV is an even number. The factor Nl/[(N/2)l] 2 is the number of 
charge-neutral terms contained in the binomial expansion of the left-hand side. The 
partition function (3.128) may thus be written as [p4|] 



* = £ 



(2z) 



N 



Nl 



N=0 

D0exp 



D^exp 



2z cos(q4>) 



N 



(3.131) 



where in the final form we recognize the sine-Gordon model. This is the dual theory we 
were seeking. Contrary to the original formulation (3.1 19), which contains the vortices 
as singular objects, the dual formulation has no singularities. To see how the vortices 
and the Kosterlitz-Thouless phase transition are represented in the dual theory we note 
that the field equation of the auxiliary field now reads 



iTV 2 (j> = 2irzq (e 



iqcji 



-iq<t> 



)• 



(3.132) 



On comparison with the previous field equation ( 3.129| ), it follows that the right-hand 
side represents the charge density of the Coulomb gas. In terms of the scalar potential 
$, Eq. (3.132) becomes the Poisson-Boltzmann equation 



= -2irq (ze 



ze 



0q$>\ 



(3.133) 



describing, at least for temperatures above the Kosterlitz-Thouless temperature, a plasma 
of positive and negative charges with density n±, 



n± 



ze 



(3.134) 



respectively. The fugacity z is the density at zero scalar potential. (It is to r ecalled that 
we suppress factors of a denoting the core size of the vortices.) Equation ( 3.133 ) is a 
self-consiste nt equ ation for the scalar potential $ giving the spatial distribution of the 
charges via ( 3.134 ). It follows from this argument that the interaction term 2z cos(q<fi) 
of the sine-Gordon model represents a plasma of vortices. 

The renormalization group applied to the sine-Gordon model reveals that at the 
Kosterlitz-Thouless temperature Tkt = \q 2 there is a phase transition between a low- 
temperature phase of tightly bound neutral pairs and a high-temperature plasma phase 
of unbound vortices [Wfl]. In the low-temperature phase, the (renormalized) fugac- 
ity scales to zero in the large-scale limit so that the interaction term, representing the 
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plasma of unbound vortices, is suppressed. The long-distance behavior of the low- 
temperature phase is therefore well described by the free theory (V</>) 2 /47r/3, repre- 
senting a gapless mode — the so-called Kosterlitz-Thouless mode. This is the superfluid 
state. The expectation value of a single vortex vanishes because in this gapless state its 
energy diverges in the infrared. 

An important characteristic of a charged plasma is that it has no gapless excita- 
tions, the photon being transmuted into a massive plasmon. To see this we assume that 
g$ << T, so that sinh(/3(7$) sa (3q&. In this approximation, the Poisson-Boltzmann 
equation (3.133) can be linearized to give 



(V 2 - m 2 })* = 0, = 4ir/3zq 2 



(3.135) 



This shows us that, in contradistinction to the low-temperature phase, in the high- 
temperature phase, the scalar potential describes a massive mode — the plasmon. In 
other words, the Kosterlitz-Thouless mode acquires an energy gap ran- Since it pro- 
vides the high-temperature phase with an infrared cutoff, isolated vortices have a finite 
energy now and accordingly a finite probability to be created. This Debeye mechanism 
of mass generation for the photon should be distinguished from the Higgs mechanism 
which operates in superconductors (see below) and which also generates a photon mass. 

Another property of a charged plasma is that it screens charges. This so-called 
Debeye screening may be illustrated by adding an external charge to the system. The 
linearized Poisson-Boltzmann equation (3.135) then becomes 



(V 2 - m£)$(x) = -27rq <5(x), 



(3.136) 



with q the external charge which we have placed at the origin. The solution of this 
equation is gi ven by $(x) = qoKo(rnjy\x\) with Kq a modified Bessel function. The 
mass term in ( 3.136 ) is {2tt times) the charge density induced by the external charge, 
i.e., 

l(x) =-— q mlK {moM). (3.137) 



/find I 



2tt- 



By integrating this density over the entire system, we see that the total induced charge 
/ x Pind = —qa completely screens the external charge — at least in the linear approx- 
imation we are using here. The inverse of the plasmon mass is the so-called Debeye 
screening length. 

To see that the sine-Gordon model gives a dual description of a superfluid film we 
cast the field equation (3.129) in the form 



iTV 2 <f> = -mqV x v s 



(3.138) 



where we employed Eq. ( |3.98[ ). On integrating this equation, we obtain up to an irrel- 
evant integration constant 



iTdi4> = -qe ij (d J Lp - ipj). 



(3.139) 



This relation, involving the antisymmetric Levi-Civita symbol, is a typical one between 
dual variables. It also nicely illustrates that although the dual variable is a regular 



©Amstex 



3.6 Dual Theory 



35 



field, it nevertheless contains the information about the vortices which in the original 
formulation are described via the singular vortex gauge field (p p . 

Given this observation it is straightforward to calculate the current-current correla- 
tion function («7,-(k)</j(— k)), with 



g = PsV s 



the mass current. We find 

( ft (k) 5j (-k)> 



2tt/3 ; 



rttktjikkh (0(k)0(-k)) , 



where the average is to be taken with respect to the partition function 



Z 



D<^exp 



1 



(3.140) 



(3.141) 



(3.142) 



which is obtained from ( p . 1 3 1| ) by setting the interaction term to zero. We obtain in 



this way the standard expression for a superfluid 

Ps 1 



( 5j (k) 5j (-k)> = 



k 2 



(3.143) 



The 1 /k 2 reflects the gaplessness of the </>-field in the low-temperature phase, while the 
combination (5y k 2 — kikj arises because the current is divergent free, V • g(x) = 0, or 
k-g(k)=0. 
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Chapter 4 

Superconductivity 



In this chapter we shall demonstrate a close connection between the Bogoliubov the- 
ory of superfluidity discussed in the previous chapter and the strong-coupling limit of 
the BCS theory of superconductivity. The phase-only effective theory governing the 
superconducting state is derived. It is also pointed out that a superconducting film at 
finite temperature undergoes a Kosterlitz-Thouless phase transition. 



4.1 BCS Theory 

Our starting point is the famous microscopic model of Bardeen, Cooper, and Schrieffer 
(BCS) defined by the Lagrangian [0 

C = V T * Wo - £(-iV)ty T + VI [ido - £(-*V)M - A € 1>l ^ 

■■= £o+4, (4-1) 

where C\ = — Ao^* ip* ipi ifif is a contact interaction term, representing the effective, 
phonon mediated, attraction between electrons with coupling constant Aq < 0, and 



£0 is the remainder. In (ft-.ip, the field t/'TU) is an anticommuting field describing the 
electrons with mass m and spin up (down); £(— iV) = e(— iV) — Mo, with e(— zV) = 
— V 2 /2m, is the kinetic energy operator with the chemical potential /j,q subtracted. 



The Lagrangian (4.1) is invariant under global U(l) transformations. Under such a 



transformation, the electron fields pick up an additional phase factor 

t/v -» e ia Va (4.2) 
with a =t, J. and a a constant. Notwithstanding its simple form, the microscopic 



model (4.1) is a good starting point to describe BCS superconductors. The reason is 
that the interaction term allows for the formation of Cooper pairs which below a crit- 
ical temperature condense. This results in a nonzero expectation value of the field A 
describing the Cooper pairs, and a spontaneous breakdown of the global U(l) symme- 
try. This in turn gives rise to the gapless Anderson-Bogoliubov mode which — after 
incorporating the electromagnetic field — lies at the root of most startling properties of 
superconductors | 
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To obtain the effective theory governing the Anderson-Bogoliubov mode, let us 
integrate out the fermionic degrees of freedom. To this end we introduce Nambu's 



notation and rewrite the Lagrangian ( \. 1 ) in terms of a two-component field 



In this notation, £ becomes 







La - v I o td + c(-tV) 1 v ' 



(4.3) 



(4.4) 



where we explicitly employed the anticommuting character of the electron fields and 
neglected terms which are a total derivative. The partition function, 



Z = J Thp^Thpexp (ij Cj , 



(4.5) 



must for our purpose be manipulated in a form bilinear in the electron fields. This is 
achieved by rewriting the quartic interaction term as a functional integral over auxiliary 
fields A and A* (for details see Ref. [El): 



exp f-i X J ip* tp* V»4. V'T 



(4.6) 



DA*DAexp 



/ / | A* ip± + if)* ip^ A — — A* A 



where, as always, an overall normalization factor is omitted. Classically, A merely 
abbreviates the product of two electron fields 



A = A ^j.^T- 



(4.7) 



It would therefore be more appropriate to give A two spin labels Aj,f. Since tp^ and Vx 
are anticommuting fields, A is antisymmetric in these indices. Physically, it describes 
the Cooper pairs of the superconducting state. 



By employing ( |4.6[ ), we can cast the partition function in the desired bilinear form: 

(4.8) 



Z = J D^DtA J DA*DA cxp (j- J A*A^j 



x exp 



t ( ido-^-iV) 
-A* 



-A 
z9 + C(-«V) 



Changing the order of integration and performing the Gaussian integral over the Grass- 
mann fields, we obtain 

Z = J DA*DA exp U,S eS [A*, A] + y J A*A^j , (4.9) 
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where S e g is the one-loop effective action which, using the identity Det(^4) = exp[Tr 
ln(A)], can be cast in the form 

S,„[A.,A] = - i I i ln(«>:i! P) M ;f (p) ), (4-10) 
where po = ida an d £(p) = £ (p) — Mo. with e(p) = p 2 /2m. 



In the mean-field approximation, the functional integral (4.9) is approximated by 
the saddle point: 

Z = exp (iS eS \A* mU A mf ] + j- J A^A m£ ^ , (4.11) 
where A m f is the solution of mean-field equation 

dA*(x) Aq 



If we assume the system to be spacetime independent so that A m f (x) = A, Eq. (4.12) 
yields the celebrated BCS gap [[47| equation: 



1 . f 1 



A " J k k%-E 2 (k)+iT] 

1 f 1 



2 7 k E(k)' 



(4.13) 



where -q is an infinitesimal positive constant that is to be set to zero at the end of the 
calculation, and 



E(k) = ^*(k) + |A|2 (4.14) 

is the spectrum of the elementary fermionic excitations. If this equation yields a solu- 
tion with A 7^ 0, the global U(l) symmetry (4.2) is spontaneously broken since 

A — > e 2iQ A ^ A (4.15) 

under this transformation. The factor 2 in the exponential function arises because A, 
describing the Cooper pairs, is built from two electron fields. It satisfies Landau's 
definition of an order parameter as its value is zero in the symmetric, disordered state 
and nonzero in the state with broken symmetry. It directly measures whether the U(l) 
symmetry is spontaneously broken. 

In the case of a spacetime-independent system, the effective action ( 4. 1C ) is readily 
evaluated. Writing 

Po-C(p) "A \ = f Po-UP) \ _/ A 

-A* p +a(p) J \ p + c(p)y v a* 

(4.16) 

and expanding the second logarithm in a Taylor series, we recognize the form 

S eff [A % A] = -*Tr m ^°-^) ° lA ' 2 ^ 



Po+fo) j \ pi-e{p))' 

(4.17) 
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up to an irrelevant constant. The integral over the loop energy ko gives for the corre- 
sponding effective Lagrangian 

C eS = J [E(k) - C(k)] . (4.18) 

To this one-loop result we have to add the tree term |A| 2 /Ao. Expanding E(k) in A, 
we see that the effective Lagrangian also contains a term quadratic in A. This term 
amounts to a renormalization of the coupling constant; we find to this order for the 
renormalized coupling constant A: 

1 1 1 ' 1 (4..9, 



A A„ 2 i t |«k)| 

where it should be remembered that the bare coupling constant Ao is negative, so that 
there is an attractive interaction between the fermions. We shall analyze this equation 
later on, for the moment it suffice to note that we can distinguish two limits. One, 
the limit where the bare coupling constant is taken to zero, Ao — > 0~, which is the 
famous weak-coupling BCS limit. Second, the limit where the bare coupling is taken 
to minus infinity Ao — > — oo. This is the strong-coupling limit, where the attractive 
interaction is such that the fermions form tightly bound pairs [|50||. These composite 
bosons have a weak repulsive interaction and can undergo Bose-Einstein condensation 
(see succeeding section). 

Since there are two unknowns contained in the theory, viz., A and fio, we need a 



second equation to determine these variables in the mean-field approximation [ 50 1 . To 



find the second equation we note that t he ave rage fermion number N, which is obtained 



by differentiating the effective action ( |4.10| ) with respect to /i 



N=^, (4.20) 

is fixed. If the system is spacetime independent, this reduces in the one-loop approxi- 
mation to 

n = -itr / Go(fc)r 3 , (4.21) 

Jk 

where n = N/V, with V the volume of the system, is the constant fermion number 
density, t% is the diagonal Pauli matrix in Nambu space, 

T3 = ( I ? ), (4.22) 



-1 



and Go(k) is the Feynman propagator, 



1 f k e lko l + £(k) A \ 

fcg - E 2 (k) + i V \ A* k e-*>" - £(k) )' 
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Here, rj is an infinitesimal positive constant that is to be set to zero at the end of the 
calculation. The exponential functions in the diagonal elements of the propagator are 
an additional convergence factor needed in nonrelativistic theories [|1J). If the integral 



over the loop energy fco in the particle number equation ( |4.21| ) is carried out, it takes 
the familiar form 



IX 1 -m) 



The two equations (4.13) and ( 4.24 ) determine A and /iq ■ They are usually evaluated in 
the weak-coupling BCS limit. However, as was first pointed out by Leggett [|50|l, they 
can also be easily solved in the strong-coupling limit (see succeeding section), where 
the fermions are tightly bound in pairs. More recently, also the crossover between 
the weak-coupling BCS limit and the strong-coupling composite boson limit has been 
studied in detail [gll g ||]. 

We are now in a position to derive the effective theory governing the Anderson- 
Bogoliubov mode. To this end we write the order parameter A m f as 

A ml (x) = Ae 2 ^W, (4.25) 



where A is a spacetime-independent solution of the mean-field equation (4.12) and 
<p{x) represents the Anderson-Bogoliubov mode, i.e., the Goldstone mode of the spon- 
taneously broken U(l) symmetry. This approximation, where the phase of the order 
parameter is allowed to vary in spacetime while the modulus is kept fixed, is called 
the London limit. This limit is relevant for our discussion of the zero-temperature 
superconductor-to-insulator phase transition in Ch. || because this transition is driven 
by phase fluctuations; the modulus of the order parameter remains finite and constant 
at the transition. The critical behavior can thus be studied with this effective theory for- 
mulated solely in terms of the phase field. We proceed by decomposing the Grassmann 
field as, cf. @] 

Mx) = e Mx) xAx) (4-26) 

and substituting the specific form ( |4.25|) of the order parameter in the partition function 
(4.8). Instead of the effective action ( |4.To| ) we now obtain 



c t , ( Po - d <p - £(p + V(f) -A \ 

o G ff = — ilrlri n . j- / rT \ ) (4.27) 

V Po + do<P + £(p - V<p) J 

where the derivative d^ip of the Goldstone field plays the role of an Abelian gauge field. 
This expression can be handled with the help of the derivative expansion outlined in 
Sec. [2~2| , to yield the phase-only effective theory. We shall not give any details here and 
merely state the result ijjjt], that the effective theory is again of the form (3.22). 



4.2 Composite Boson Limit 

In this section we shall investigate the strong-coupling limit of the pairing theory. In 
this limit, the attractive interaction between the fermions is such that they form tightly 
bound pairs of mass 2m. To explicate this limit in arbitrary dimension d, we swap the 
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bare coupling constant for a more convenient parameter — the binding energy e a of a 
fermion pair in vacuum Both parameters characterize the strength of the contact 
interaction. To see the connection between the two, let us consider the Schrodinger 
equation for the problem at hand. In reduced coordinates it reads 



V 2 

+ A S x) 

to 



V(x) = 



(4.28) 



where the reduced mass is to/2 and the delta-function potential, with Ao < 0, repre- 
sents the attractive contact interaction d in (4T). We stress that this is a two-particle 
problem in vacuum; it is not the famous Cooper problem of two interacting fermions 
on top of a filled Fermi sea. The equation is most easily solved by Fourier transforming 
it. This yields the bound-state equation 



V(k) 



A„ 



k 2 /m 



-^(0), 



or 



Ao Jk 



(4.29) 



(4.30) 



/ k k 2 /TO + e a ' 

This equation allows us to replace th e cou pling constant with the binding energy e a . 
When substituted in the gap equation (4.13), the latter becomes 



1 



1 



, k k 2 /TO + e a 2j k E(k) ■ 
By inspection, it is easily seen that this equation has a solution given by [p0|] 



0, Mo 



(4.31) 



(4.32) 



where it should be noted that the chemical potential is negative here. This is the strong- 
coupling limit. To appreciate the physical significance of the specific value found for 
the chemical potential in this limit, we note that the spectrum (q) of the two-fermion 
bound state measured relative to the pair chemical potential 2^i reads 



£b(q) 



q 

4m 



(4.33) 



The negative value for no found in ( |4.32 ) is precisely the condition for a Bose-Einstein 
condensation of the composite bosons in the q — state. 

To investigate this limit further, we consider the effective action ( 1.1 C| ) and expand 
A(x) around a constant value A satisfying the gap equation (4.13), 



AO) =A + A{x). 



(4.34) 



We obtain in this way, 



i Tr 



^ I 

i=i 



G (P) 





A* 



A 





(4.35) 
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where Go is given in (4.23). We are interested in terms quadratic in A. Employing the 

(4.36) 



derivative expansion outlined in Sec. 2.2, we find 

1 



1 



x {A 2 A* A* + [po + e(p)][po + 9o - e(p - q)]AA* 
+A* 2 AA + [p - f(p)][po + qo + f(p - q)]A*A} ; 

where g M = i^,. It is to be recalled here that the derivative p^ operates on everything 
to its right, while operates only on the first object to its right. Let us for a moment 
ignore the derivatives in this expression. After carrying out the integral over the loop 
energy ko and using the gap equation ( \. 1 3 ), we then obtain 

1 



£( 2 >(0) 



£ 3 (k) 



A 2 A* 2 



A* 2 A 2 + 2|A| 2 |A| 2 



(4.37) 



In the composite boson limit A — * 0, so that the spectrum (1 14) of the elementary 
fermionic excitations can be approximated by 



E(k) a e(k) 



The remaining integrals in ( 4.37 ) then become elementary, 

1 = 4T(3-d/2) d/2 rf/2 _ 3 
E 3 (k) (47r) d / 2 



(4.38) 



(4.39) 



We next consider the terms involving derivatives in (4.36). Following Ref. [ p2[ | we 
set A to zero here. The integral over the loop energy is easily carried out, with the 
result 



k qo — k 2 /m + 2^ — q 2 /4m 
1 



-qo — k 2 jm + 2/io — q 2 /4m 



AA* 
A*A. 



(4.40) 



The integral over the loop momentum k gives in the strong-coupling limit using di- 
mensional regularization 



qo — k 2 /m — e a — q 2 /4m 
or expanded in derivatives 



(47r) d / 2 V 



qo 



(4.41) 



1 



k qo - k 2 /m - e a - q 2 /4m 

T(l - d/2) d/2 d/2 _ 1 _ T(2 - d/2) d /2 d/2-2 
( 47r )d/2 e - ( 47r )d/2 e « 



(4.42) 



90 



q 

4m 
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The first term at the right-hand side yields as contribution to the effective theory 



(2) _ r(l-rf/2) d/2 rf/2-l|X|2 

A ~~ (47r) d / 2 a 11 ' 



(4.43) 



To this we have to add the contribution | A| 2 /Ao coming from the tree potential, i.e., the 
last term in the partition function (4.9). But this combination is no other than the one 
needed to defined the renormalized coupling constant via (4.19), which in the strong- 
coupling limit reads explicitly 



I _ i m-d/2) d/ 2 d/2 _ x 

A A (4^) d / 2 a ' 



(4.44) 



In other words, the contribution ( 1.43| ) can be combi ned w ith the tree contribution to 
yield the term |A| 2 /A. E xpanding the square root in ( 4.41 ) in powers of the derivative 
using the value ( 4.32 ) for the chemical potential, and pasting the pieces together, we 
obtain for the terms quadratic in A [p2Q, 



lT(2-d/2) 



d/2 £ d/2-2| r t M (g)#, # = 



A 

A* 



2 (47T) d / 2 

where M (q) is the 2x2 matrix, 

M (q) = 

q - q 2 /4m - (2 - d/2)\A\ 2 /e a -(2 - d/2)A 2 /e a 

-(2 - d/2)A* 2 /e a -q - q 2 /4m - (2 - d/2)\A\ 2 /e a 



(4.45) 



(4.46) 



This Lagrangian is precisely of th e form found in ( |3.8| ) describing an interacting Bose 
gas. On comparing with Eq. (3.9 ), we conclude that the composite bosons have — as 
expected — a mass = 2m twice the fermion mass m, and a small chemical potential 



M0l b = (2-d/2)i 



(4.47) 



From ( 4.46| ) one easily extracts the Bogoliubov spectrum and the velocity Co of the 
sound mode it describes, 



m h me a 



Mo.b 



Also the number density no.b of condensed composite bosons, 

^o, b = ^|^^ /2 - 2 |A| 2 
as well as the weak repulsive interaction Ao.b between the composite bosons, 



A ,b = (4 7 r) d / 1 



l-d/4 e 1 ^ 2 
r(2-d/2) m d / 2 



(4.48) 



(4.49) 



(4.50) 
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follow immediately. We in this way have explicitly demonstrated that the BCS theory 
in the composite boson limit maps onto the Bogoliubov theory. 

In concluding this section, we remark that in d = 2 various integrals we encoun- 
tered become elementary for arbitrary values of A. For example, the gap equation 



(4.31 ) reads explicitly in d = 2 



Hi + |A| 2 - m , (4.51) 



while the particle number equation (4.24) becomes 

rii 



Since in two dimensions, 



«=^(vVo + |A| 2 + Mo)- (4.52) 



7T- = — cf, (4-53) 



2n 7T 



with fcp and ef = k\j1m the Fermi momentum and energy, the two equations can be 



combined to yield [ 57 1 



tg _ 2 y//^ + |A| 2 -Mo (454) 
e F VMo + IAP + Zio 

The composite boson limit we have been discussing in this section is easily retrieved 
from these more general equations. Also note that in this limit, h — 2no,t» while the 
renormalization of the coupling constant takes the same form as for an interacting Bose 
gas 

1 1 f 1 m \ „ «x 



cf. (3.32) 



4.3 Dual Theory 

We now turn to the dual description of a superconducting film at finite temperature. 



We thereto minimally couple the model of Sec. 3.5 to a magnetic field described by 
the magnetic vector potential A. For the time being we ignore vortices by setting the 
vortex gauge field (p p to zero. The partition function of the system then reads 

Z = J Dtp J DA 3(A) exp f-/3 j Tij , (4.56) 

where 5(A) is a gauge-fixing factor for the gauge field A, and TC is the Hamiltonian 

H = ipsv 2 |i(Vx A) 2 (4.57) 

with 

v s = — (S7ip-2eA). (4.58) 
m 
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The double charge 2e stands for the charge of the Cooper pairs which are formed at the 
bulk transition temperature. The functional integral over ip in ( |4.56 ) is easily carried 
out with the result 



Z = /DA S (A)exp{-f/_ 



(V x A) 



m 2 ^ 



d t d. 



A, 



(4.59) 



where the last term, with m 2 A — 4e 2 p s /m 2 , is a gauge-invariant, albeit nonlocal mass 
term for the gauge field generated by the Higgs mechanism. The number of degrees of 
freedom does not change in the process. This can be seen by noting that a gapless gauge 
field in two dimensions represents no physical degrees of freedom. (In Minkowski 
spacetime, this is easily understood by recognizing that in 1 + 1 dimensions there is 
no transverse direction.) Before the Higgs mechanism took place, the system therefore 
contains only a single physical degree of freedom described by ip. This equals the 
number of degrees of freedom contained in (4.59). 



We next introduce an auxiliary field h to linearize the first term in (4.59) 



exp 



-§/<VxA)= 



= / Dhexp 



— I h 2 



h(V x A) 



(4.60) 



and integrate out the gauge-field fluctuations [with a gauge-fixing term (l/2a)(V ■ 
A) 2 ]. The result is a manifestly gauge-invariant expression for the partition function in 
terms of a massive scalar field h, representing the single degree of freedom contained 
in the theory: 



Z = Dh exp 



1 



1 



(4.61) 



To understand the physical significance of this field, we note from (4.60) that it satisfies 
the field equation 

h = if3V x A. (4.62) 

That is, the fluctuating field h represents the local magnetic induction, which is a scalar 
in two space dimensions. Equation (4.61) shows that the magnetic field has a finite 



penetration depth Al = I/vtla- In contrast to the original description wh ere th e func- 
tional integral runs over the gauge potential, the integration variable in ( 4.61 ) is the 
physical field. 

We next include vortices. The penetration depth Al provides the system with an 
infrared cutoff so that a single magnetic vortex in the charged theory has a finite energy. 
Vortices can therefore be thermally activated. This is different from the superfluid 
phase of the neutral model, where the absence of an infrared cutoff permits only tightly 
bound vortex-antivortex pairs to exist. We expect, accordingly, the superconducting 
phase to describe a plasma of vortices, each carrying one magnetic flux quantum ±7r/e. 
The partition function now reads 



00 r Af++JV_ 

y - — 



Hj J Dp J DAS(A) exp ( 4 -63) 
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where z is the fugacity, i.e., the Boltzmann factor associated with the vortex core en- 
ergy. The velocity appearing in the Hamiltonian (4.57) now includes the vortex gauge 
field 



v s = — (V<p-2eA- <p p ). 



(4.64) 



This field can be shifted from the first to the second term in the Hamiltonian (4.57) by 
applying the transformation A — > A — <p p /2e. This results in the shift 



VxA^VxA-B b 



with the plastic field 



IV = -$ ^w Q <5(x-x") 



(4.65) 



(4.66) 



representing the magnetic flux density. Here, $o = 7r/e is the elementary fl ux qu an- 
tum. Repeating the steps of the previous paragraph we now obtain instead of ( 4.61 ) 



Z 



^ N+\N-\ 

N ± =0 ^ 



n 



Dh exp 



1 



1 



(Vhf 



+ i 



B p h' { 
(4.67) 



where h represents the physical local magnetic induction h 

h = if3(V x A-B p ) =i/3h. 



(4.68) 



The field equation for h obtained from (4.67 ) yields for the magnetic induction: 

- V 2 h + m\h = m 2 A B p : 



(4.69) 



which is the familiar equation in the presence of magnetic vortices. 

The last term in (4.67) shows that the charge g with which a magnetic vortex cou- 
ples to the fluctuating /i-field is the product of an elementary flux quantum (contained 
in the definition of B F ) and the inverse penetration depth ttia = 1/ Al, 



g = $ m A . 



(4.70) 



For small fugacities the summation indices N + and 7V_ can be restricted to the values 
0, 1 and we arrive at the partition function of the massive sine-Gordon model [p8[| 



Z= / D/m-. m ,(- / {-L 



rr> 



- 2ZCOS ( $n/l 



(4.71) 



This is the dual formulation of a two-dimensional superconductor. The magne tic vor- 
tices of unit winding number w a = ±1 turned the otherwise free theory ( 4.61 ) into an 
interacting one. 

The final form ( 4.71 ) demonstrates the rationales for going over to a dual theory. 
First, it is a formulation directly in terms of a physical field representing the local 
magnetic induction. There is no redundancy in this description and therefore no gauge 
invariance. Second, the magnetic vortices are accounted for in a nonsingular fashion. 
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This is different from the original formulation of the two-dimensional superconductor 
where the local magnetic induction is the curl of an unphysical gauge potential A, and 
where the magnetic vortices appear as singular objects. 

Up to this point we have discussed a genuine two-dimensional superconductor. As 
a model to describe superconducting films this is, however, not adequate. The reason 
is that the magnetic interaction between the vortices takes place mostly not through 
the film but through free space surrounding the film where the photon is gapless. This 
situation is markedly different from a superfiuid film. The interaction between the 
vortices there is mediated by the Kosterlitz-Thouless mode which is confined to the 
film. A genuine two-dimensional theory therefore gives a satisfactory description of a 
superfiuid film. 

To account for the fact that the magnetic induc tion is not confined to the film and 
can roam in outer space, the field equation ( |4.69 ) is modified in the following way 



- V 2 /i(x ± , X3 ) + -^-5 d (x 3 )h(x ± ,x 3 ) = -^-S d (x 3 )B p ( X ). (4.72) 

Here, 1/Aj_ = dm 2 A , with d denoting the thickness of the superconducting film, is 
an inverse length scale, xi denotes the coordinates in the plane, h the component of 
the induction field perpendicular to the film, and S d (x 3 ) is a smeared delta function of 
thickness d along the X3-axis 

P / n f = for b 3 | > d/2 

^0 for 1^1 < d/2 ' (473) 



The reason for including the smeared delta function at the right-hand side of ( 4.72 ) is 
that the vortices are confined to the film. The delta function in the second term at the 
left-hand side is included because this term is generated by screening currents which 
are also confined to the film. 

To be definite, we consider a single magnetic vortex located at the origin. The 



induction field found from (4.72) reads 



fc(xi '° )= £r^T+fc^ j ° (ff|xi|) ' (4 - ?4) 

with Jo the Oth Bessel function of the first kind. At small distances from the vortex 
core (Xj_q » 1) 

ft(xx,0)~ $ ° (4.75) 
47TA Jxjl 



while far away (X±q << 1) 



7T XJ_ 



This last equation shows that the field does not exponentially decay as would be the 
case in a genuine two-dimensional system. The reason for the long range is that most 
of the magnetic interaction takes place in free space outside the film where the photon 
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is gapless. If, as is often the case, the length Aj_ = l/dm 2 A is much larger than the 
sample size, it can be effectively set to infinity. In this limit, the effect of the mag- 
netic interaction diminishes, as can be seen from ( |4.75| ), and the vortices behave as 
in a superfluid film. One therefore expects a superconducting film to also undergo a 
Kosterlitz-Thouless transition at a temperature Tkt characterized by an unbinding of 
vortex-antivortex pairs. The first experiment to study this possibility was carried out in 
Ref. [[HJ. Because the transition temperature Tkt is well below the bulk temperature 
T c where the Cooper pairs form, the energy gap of the fermions remains finite at the 
critical point [|2|| . This prediction has been corroborated by experiments performed by 
Hebard and Paalanen on superconducting films [|63[]. For temperatures Tkt < T < T c , 
there is a plasma of magnetic vortices which disorder the superconducting state. At 
Tkt vortices and antivortices bind into pairs and algebraic long-range order sets in. 
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Chapter 5 

Fractional Quantized Hall 
Effect 



The nonrelativistic |</>| 4 -theory describing an interacting Bose gas is also of importance 
for the description of the fractional quantized Hall effect (FQHE). As a function of 
the applied magnetic field, this two-dimensional system undergoes a zero-temperature 
transition between a so-called quantum Hall liquid, where the Hall conductance is 
quantized in odd fractions of e 2 /27r, or, reinstalling Planck's constant, e 2 /h, and an 
insulating phase. Here, the nonrelativistic |</>| 4 -fheory describes — after coupling to a 
Chern-Simons term — the original electrons bound to an odd number of flux quanta. 
The Hall liquid is in this picture characterized by a condensate of composite particles. 

5.1 Chern-Simons-Ginzburg-Landau Theory 

The fractional quantized Hall effect (FQHE) is the hallmark of a new, intrinsically 
two-dimensional condensed-matter state — the quantum Hall liquid. Many aspects of 
this state are well understood in the framework of the quantum-mechanical picture 



developed by Laughlin [64|. Considerable effort has nevertheless been invested in 
formulating an effective field theory which captures the essential low-energy, small- 
momentum features of the liquid. A similar approach in the context of superconduc- 
tors has proven most successful. Initially, only the phenomenological model proposed 
by Ginzburg and Landau [ |65| | in 1950 was known here. Most of the fundamental 
properties of the superconducting state such as superconductivity — the property that 
gave this condensed-matter state its name, Meissner effect, magnetic flux quantization, 
Abrikosov flux lattice, and Josephson effect, can be explained by the model. The mi- 
croscopic theory was given almost a decade later by Bardeen, Cooper, and Schrieffer 
[p7[]. Shortly here after, Gorkov [|6[ made the connection between the two approaches 
by deriving the Ginzburg-Landau model from the microscopic BCS theory, thus giving 
the phenomenological model the status of an effective field theory. 

A first step towards an effective field theory of the quantum Hall liquid was taken 
by Girvin and MacDonald [p7[| and has been developed further by Zhang, Hansson 
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and Kivelson [pSQ, who also gave an explicit construction starting from a microscopic 
Hamiltonian. Their formulation (for a review see Ref. [|69[]) incorporates time depen- 
dence which is important for the study of quantum phase transitions. 

An essential ingredient for obtaining an effective theory of the FQHE was the iden- 
tification by Girvin and MacDonald [|7| of a bosonic operator <j) exhibiting (algebraic) 
off-diagonal long-range order of a type known to exist in two-dimensional bosonic su- 
perfiuids. They argued that this field should be viewed as an order parameter in terms of 
which the effective field theory should be formulated. To account for the incompress- 
ibility of the quantum Hall liquid they suggested to minimally couple cj> to a so-called 
statistical gauge field (ao, a) governed solely by a Chern-Simons term 

£ cs = I e 2 6»<9 a x a - e 2 6»a V x a, (5.1) 

with V x a the statistical magnetic field and 9 a constant. As we will see below, the 
gapless Bogoliubov spectrum of the neutral system changes as a result of this coupling 
into one with an energy gap ^j8|], thus rendering the charged system incompressible. 

Because of the absence of a kinetic term (the usual Maxwell term), the statistical 
gauge field does not represent a physical degree of freedom. In a relativistic setting, a 
Maxwell term is usually generated by quantum corrections so that the statistical gauge 
field becomes dynamical at the quantum level. The quantum theory then differs quali- 
tatively from the classical theory. On the other hand, as we shall see below, this need 
not be the case in a nonrelativistic setting. That is to say, the Ansatz of the absence of 
a Maxwell term is here not necessarily obstructed by quantum corrections. 

The effective theory of the quantum Hall liquid is given by the so-called Chern- 
Simons-Ginzburg-Landau (CSGL) Lagrangian [p8| 



C = ipDot - ^-\m 2 + Mol^l 2 - Ao|0| 4 + Ccs- (5.2) 

The covariant derivatives Dq = do + ieAo + ieao and D = V — ieA — iea give a 
minimal coupling to the applied magnetic and electric field described by the gauge field 
(A , A) and also to the statistical gauge field. For definiteness we will assume that our 
two-dimensional sample is perpendicular to the applied magnetic field, defining the z- 
direction, and we choose the electric field to point in the x-direction. The charged field 
4> represents the Girvin-MacDonald order parameter describing the original electrons 
bound to an odd number 21 + 1 of flux quanta. To see that it indeed does, let us consider 
the field equation for ao: 

|0| 2 = _ e (9V x a. (5.3) 
The simplest solution of the CSGL Lagrangian is the uniform mean-field solution 

\4>\ 2 =n, a = -A, a = -Ao = 0, (5.4) 

where n indicates the constant fermion number density. The statistical gauge field is 



seen to precisely cancel the applied field. The constraint equation (5.3) then becomes 



e6H, (5.5) 
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with H the applied magnetic field. Now, if we choose O^ 1 = 2ir(2l + 1), it follows on 
integrating this equation that, as required, with every electron there is associated 21 + 1 
flux quanta: 

N =2iTT N ®> (5 - 6) 

where N® = $/ $0, w ith $ = J x H the magnetic flux, indicates the number of flux 



quanta. Equation (5.5) implies an odd-denominator filling factor vh which is defined 
by 

n 1 

"* = I/*; = a+T (5 ' 7) 



The coupling constant Ao (> 0) in ( p.2| ) is the strength of the repulsive contact 
interaction between the composite particles, and fiQ is a chemical potential introduced 
to account for a finite number density of composite particles. 

It is well known from anyon physics that the inclusion of the Chern-Simons term 
changes the statistics of the field (f> to which the statistical gauge field is coupled [|70[|. 
If one composite particle circles another, it picks up an additional Aharonov-Bohm fac- 
tor, representing the change in statistics. The binding of an odd number of flux quanta 
changes the fermionic character of the electrons into a bosonic one for the composite 
particles, allowing them to Bose condense. The algebraic off-diagonal long-range or- 
der of a quantum Hall liquid can in this picture be understood as resulting from this 
condensation. Conversely, a flux quantum carries 1/(21 + l)th of an electron's charge 
Q, and also 1/(21 + l)th of an electron's statistics 

The defining phenomenological properties of a quantum Hall liquid are easily 
shown to be described by the CSGL theory [ |58"i p$\ . From the lowest-order expres- 
sion for the induced electromagnetic current one finds 

ej, = ^- = = ^ = < - d jao ) = eH^E h (5.8) 



with E the applied electric field and where we have written the Lagrangian (5.2) as 
a sum C — Ccf, + Ccs- It follows that the Hall conductance a xy is quantized in odd 
fractions of e 2 /2ir, or, reinstalling Planck's constant, e /ft. This result can also be 
understand in an intuitive way as follows. Since the composite particles carry a charge 
e, the applied electric field gives rise to an electric current 

dN 

I=e— (5.9) 
at 

in the direction of E, i.e., the ir-direction. This is not the end of the story because the 
composite objects carry in addition to electric charge also 21 + 1 flux quanta. When the 
Goldstone field ip encircles 21 + 1 flux quanta, it picks up a factor 2ir for each of them 

<j> V-(/? = 2tt(2Z + 1). (5.10) 

Now, consider two points across the sample from each other. Let the phase of these 
points initially be equal. As a composite particle moves downstream, and crosses the 
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line connecting the two points, the relative phase Aip between them changes by 2ir(2l + 
1). This phase slippage [|l5| leads to a voltage drop across the sample given by 

Vk = -d Ap = (2Z + 1)$ ^-, (5.11) 
e dr 

where the first equation can be understood by recalling that due to minimal coupling 
doip — + dotp + eAo. For the Hall resistance we thus obtain the expected value 

Vh 2tt 
p xy = ^ = (2l + l) — . (5.12) 

If the CSGL theory is to describe an incompressible liquid, the spectrum of the 
single-particle excitations must have a gap. Without the coupling to the statistical gauge 



field, the spectrum is given by the gapless Bogoliubov spectrum (3.13). To obtain the 



single-particle spectrum of the coupled theory, we integrate out the statistical gauge 



field. The integration over ao was shown to yield the constraint (5.3) which in the 
Coulomb gauge V • a = is solved by 

a i = ^e ij ^\4>\ 2 . (5.13) 
The integration over the remaining components of the statistical gauge field is now 



simply performed by substituting (5.13) back into the Lagrangian. The only nonzero 
contribution arises from the term — e 2 |0| 2 a 2 /2m. The spectrum of the charged system 
acquires as a result an energy gap ui c 



E(k) = v / ^ c 2 + e 2 (k) + 2 Mo e(k), (5.14) 

with cu c = /io/2#?7jAo. To lowest order, the gap equals the cyclotron frequency of a 
free charge e in a magnetic field H 

. c = f - = (5.15) 

am m 

The presence of this energy gap results in dissipationless flow with a xx = 0. 

These facts show that the CSGL theory captures the essentials of a quantum Hall 
liquid. Given this success, it is tempting to investigate if the theory can also be em- 
ployed to describe the field-induced Hall-liquid-to-insulator transitions. This will be 
done in Sec. [5j| . It should however be borne in mind that both the l/|x|-Coulomb po- 
tential as well as impurities should be incorporated into the theory in order to obtain a 
realistic description of the FQHE. The repulsive Coulomb potential is believed to play 
a decisive role in the formation of the the composite particles, while the impurities are 
responsible for the width of the Hall plateaus. As the magnetic field moves away from 
the magic filling factor, magnetic vortices will materialize in the system to make up 
the difference between the applied field and the magic field value. In the presence of 
impurities, these defects get pinned and do not contribute to the resistivities, so that 
both a xx and a xy are unchanged. Only if the difference becomes too large, the system 
reverts to an other quantum Hall state with a different filling factor. 
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This chapter is devoted to continuous phase transitions at the absolute zero of temper- 
ature; so-called quantum phase transitions. Unlike in classical phase transitions taking 
place at finite temperature and in equilibrium, time plays an important role in quantum 
phase transitions. Put differently, whereas the critical behavior of classical 2nd-order 
phase transitions is governed by thermal fluctuations, that of 2nd-order quantum tran- 
sitions is controlled by quantum fluctuations. These transitions, which have attracted 
much attention in recent years (for an introductory review, see Ref. [[jj), are triggered 
by varying not the temperature, but some other parameter in the system, like the ap- 
plied magnetic field, the charge carrier density, or the disorder strength. The quantum 
phase transitions we will be discussing here are all dominated by phase fluctuations. 



6.1 Scaling 

The natural language to describe quantum phase transitions is quantum field theory. 
In addition to a diverging correlation length £, quantum phase transitions also have a 
diverging correlation time £ 4 . They indicate, respectively, the distance and time period 
over which the order parameter characterizing the transition fluctuates coherently. The 
way the diverging correlation time relates to the diverging correlation length, 

(6.1) 

defines the so-called dynamic exponent z. It is a measure for the asymmetry between 
the time and space directions and tells us how long it takes for information to prop- 
agate across a distance £. The traditional scaling theory of classical 2nd-order phase 
transitions, first put forward by Widom [|72|], is easily extended to include the time 



dimension [|29|] because relation ( 5T) implies the presence of only one independent di- 
verging scale. Let S = K — K c , with K the parameter that drives the phase transition, 
measure the distance from the critical coupling K c . A physical observable at the abso- 
lute zero of temperature O(ko, |k|, K) can in the critical region close to the transition 
be written as 

O(k ,\klK) = i d oO(M u \k\i), (T = 0), (6.2) 
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where do is the dimension of the observable O. The right-hand side does not depend 
explicitly on K; only implicitly through £ and £ t . The closer one approaches the critical 
coupling K c , the larger the correlation length and time become. 

Since a physical system is always at some finite temperature, we have to investigate 



how the scaling law ( 5.2) changes when the temperature becomes nonzero. The easiest 
way to include temperature in a quantum field theory is to go over to imaginary time 
r = it, with r restricted to the interval < r < (3. The temporal dimension becomes 
thus of finite extend. The critical behavior of a phase transition at finite temperature 
is still controlled by the quantum critical point provided £ t < (3. If this condition 
is fulfilled, the system does not see the finite extend of the time dimension. This is 
what makes quantum phase transitions experimentally accessible. Instead of the zero- 



temperature scaling (5.2), we now have the finite-size scaling 

O(k Q ,\k\,K,T)^/3 do ^O(k /3,\k\/3 1 / z ,f3/^), (T + 0). (6.3) 
The distance to the quantum critical point is measured by the ratio /3/£ t ~ \5\ zu /T. 



6.2 Repulsively Interacting Bosons 

The first quantum phase transition we wish to investigate is the superfluid-to-Mott- 



insulating transition of interacting bosons in the absence of impurities [73 1. The transi- 
tion is described by the nonrelativistic |0| 4 -theory (3.1), which becomes critical at the 



absolute zero of temperature at some (positive) value /i c of the renormalized chemical 
potential. The Mott insulating phase is destroyed and makes place for the superfluid 
phase as /i increases. Whereas in the superfluid phase the single-particle (Bogoliubov) 
spectrum is gapless and the system compressible, the single-particle spectrum of the 
insulating phase has an energy gap and the compressibility k vanishes here. 

The nature of the insulating phase can best be understood by putting the theory on 
a lattice. The lattice model is defined by the Hamiltonian 

Hn = -t^2(ata j+1 + h.c.) + ^](-/i L n 3 + Un 2 ), (6.4) 

3 3 

where the sum . is over all lattice sites. The operator a}- creates a boson at site j 

and fij = ajctj is the particle number operator at that site; t is the hopping parameter, 
U the interparticle repulsion, and ^ is the chemical potential on the lattice. The zero- 
temperature phase diagram is as follows fl73|], In the limit t/U — > 0, each site is 
occupied by an integer number n of bosons which minimizes the on-site energy (see 



Fig. 5.1) 



e(n) = + Un 2 . (6.5) 

It follows that within the interval 2n— 1 < /.il/U < 2n+l, each site is occupied by ex- 
actly n bosons. When the chemical potential is negative, n = 0. The intervals become 
smaller when t/U increases. Within such an interval, where the particles are pinned to 
the lattice sites, the single-particle spectrum has an energy gap, and the system is in the 
insulating phase with zero compressibility, k = nT 2 dn / d [i^ = 0. Outside these inter- 
vals, the particles delocalize and can hop through the lattice. Being at zero temperature, 
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Figure 6.1: Schematic representation of the average number n of particles per site as 
function of the chemical potential fi^ at some finite value of the hopping parameter 

t < t c . 



the delocalized bosons condense in a superfiuid state. The single-particle spectrum is 
gapless here and the system compressible (k ^ 0). 

Ast/U increases, the gap in the single-particle spectrum as well as the width of the 
intervals decrease and eventually vanish at some critical value t c . For values t > t c of 
the hopping parameter, the superfiuid phase is the only phase present (see Fig. 6.2). 

The continuum model (3.1), with fi > fi c describes the condensed delocalized 
lattice bosons which are present when the density deviates from integer values (see 
Fig. 6.1 ). In the limit fi — > fi c from above, the number of delocalized bosons decreases 
and eventually becomes zero at the phase boundary fi = fi c between the superfiuid and 
insulating phases. 

Various quantum phase transitions belong to the universality class defined by the 
zero-density transition of repulsively interacting bosons. For example, itinerant quan- 
tum antiferromagnets j74| , [75] ffdj as well as lower-dimensi onal (clean) superconduc- 
tors belong to this universality class. As we have seen in Sec. 4.2, Cooper pairs become 
tightly bound composite particles in the strong-coupling limit, which are described by 
the nonrelativistic |0| 4 -theory with a weak repulsive interaction. For fi > fi c , the field 
cj) now describes the condensed delocalized Cooper pairs. When the chemical poten- 
tial decreases, the condensate diminishes, and the system again becomes insulating for 
fi < fi c [|62|. By continuity, we expect also the superconductor-to-insulator transition 
of a (clean) weakly interacting BCS superconductor to be in this universality class. The 
restriction to lower dimensions is necessary for two different reasons. First, o nly f or 
d < 2 the penetration depth is sufficiently large [see, for example, below Eq. ( 4.72 )], 
so that it is appropriate to work in the limit Al — > oo with no fluctuating gauge field 
P%. Second, in lower dimensions, the energy gap which the fermionic excitations face 
remains finite at the critical point, so that it is appropriate to ignore these degrees of 
freedom. Moreover, since also the coherence length remains finite at the critical point, 
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superfluid 



Figure 6.2: Schematic representation of the phase diagram of the lattice model (6.4) at 
the absolute zero of temperature [[73]]. 



the Cooper pairs look like point particles on the scale of the diverging correlation length 
associated with the phase fluctuations, even in the weak-coupling limit [^j. 

In the preceding chapter, we argued that the nonrelativistic |(/>| 4 -theory is also of 
importance for the description of the fractional quantized Hall effect (FQHE), where 
it describes — after coupling to the Chern-Simons term — the original electrons bound 
to an odd number of flux quanta. As function of the applied magnetic field, this two- 
dimensional system undergoes a zero-temperature transition between a quantum Hall 
liquid, where the Hall conductance is quantized in odd fractions of e 2 /2ir, and an 
insulating phase. The Hall liquid corresponds to the phase with fi > fi c , while the 
other phase again describes the insulating phase. 

It should be noted however that in most of the applications of the nonrelativistic 
1 01 4 -theory mentioned here, impurities play an important role; this will be the main 



subject of Sec. 5.4 



The critical properties of the zero-density transition of the nonrelativistic |</>| 4 - 
theory were first studied by Uzunov p^. To facilitate the discussion let us make 
use of the fact that in nonrelativistic theories the mass is — as far as critical phenomena 
concerned — an irrelevant parameter which can be transformed away. This transforma- 
tion changes, however, the scaling dimensions of the 0-field and the coupling constant 
which is of relevance to the renormalization-group theory. The engineering dimensions 
become 

[x] = -l, [t] = -2, [Mo]=2, [Ao] = 2-d, M = |d, (6.6) 

with d denoting the number of space dimensions. In two space dimensions the cou- 
pling constant Ao is dimensionless, showing that the |0| 4 -term is a marginal operator, 
and d c = 2 the upper critical space dimension. Uzunov showed that below the upper 
critical dimension there appears a non-Gaussian infrared-stable (IR) fixed point. He 
computed the corresponding critical exponents to all orders in perturbation theory and 
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Figure 6.3: A closed oriented loop. 

X+ XX+XXK 
+ XXXK + 

Figure 6.4: Ring diagrams renormalizing the vertex function of the neutral |</>| 4 -theory. 



showed them to have Gaussian values, v = |, z = 2, 77 = 0. Here, 1/ character- 
izes the divergence of the correlation length, z is the dynamic exponent, and 77 is the 
correlation-function exponent which determines the anomalous dimension of the field 
<j>. The unexpected conclusion that a non-Gaussian fixed point has nevertheless Gaus- 
sian exponents is rooted in the analytic structure of the nonrelativistic propagator at 
zero bare chemical potential (p = 0): 



> -°w- iTTg+g <6 ' 7) 

where, as before, 77 is a small positive constant that has to be taken to zero after the 
loop integrations over the energies have been carried out. The factor exp(ifco?7) is an 
additional convergence factor typical for nonrelativistic theories, which is needed for 
Feynman diagrams involving only one (/(-propagator. The rule k — > k + irj in ( |6.7| ) 
expresses the fact that in this nonrelativistic theory particles propagate only forward in 
time. In diagrams involving loops with more than one propagator, the integrals over 
the loop energy are convergent and can be evaluated by contour integration with the 
contour closed in either the upper or the lower half plane. If a diagram contains a loop 
which has all its poles in the same half plane, it consequently vanishes. Pictorially, 
such a loop has all its arrows, representing the Green functions contained in the loop, 



oriented in a clockwise or anticlockwise direction | ]78| ] (see Fig. 63). We will refer 
to them as closed oriented loops. Owing to this property most diagrams are zero. 
In particular, all self-energy diagrams vanish. The only surviving ones are the so- 



called ring diagrams which renormalize the vertex (see Fig. 6.4). Because this class of 
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diagrams constitute a geometric series, the one-loop result is already exact. The vertex 
renormalization leads to a non-Gaussian fixed point in d < 2, while the vanishing 
of all the self-energy diagrams asserts that the exponents characterizing the transition 
are not affected by quantum fluctuations and retain their Gaussian values [E4j|. These 
results have been confirmed by numerical simulations in d = 1 [ f7^ ] and also by general 
scaling arguments | 



0, 73| 



To understand the scaling arguments, let us consider the two terms in the effective 



theory (3.22) quadratic in the Goldstone field tp with m effectively set to 1 [fc 



1^2 



K,{d ip) 2 



(6.8) 



We have written this in the most general form. The coefficient p s is the superfluid mass 
density which in the presence of, for example, impurities does not equal mn — even at 
the absolute zero of temperature. The other coefficient, 



n 2 K 



— = lim n 00 (o,k), 

op k^a 



(6.9) 



with n 00 the (0 0)-component of the full polarization tensor ( 3.19 ), involves the full 
compressibility and particle number de nsity of the system at rest. This is because 
the chemical potential is according to ( 3.23 ) represented in the effective theory by 
p = —dop and 

d 2 C 

(6.10) 



dp 2 



n 2 K. 



Equation ( p.8[ ) leads to the general expression of the sound velocity 



Pa 



(6.11) 



at the absolute zero of temperature. 

Let S oc p — p c denote the distance from the phase transition, so that £ ~ \S\~ V . 
Now, on the one hand, the singular part of the free energy density / S mg arises from the 
low-energy, long-wavelength fluctuations of the Goldstone field. (Here, we adopted 
the common practice of using the symbol / for the density ft/V and of referring to it 
as the free energy density.) The ensemble averages give 



((V^) : 



((W 2 >~c 3 ~r 



(6.12) 



On the other hand, dimensional analysis shows that the singular part of the free energy 
density scales near the transition as 



fsi 



-(d+z) 



(6.13) 



Combining these hyperscaling arguments, we arrive at the following conclusions: 

ps ^ ^-(d+z-2) ^ -2 K _ £-(d-*) _ \5\(d-*>, ( 6 .14) 



The first conclusion is consistent with the universal jump (3.125) predicted by Nelson 
and Kosterlitz [45 1 which corresponds to taking z = and d = 2. Since £ ~ |<5| _i/ , 
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/sing can also be directly differentiated with respect to the chemical potential to yield 
for the the singular part of the compressibility 

2 K„: |£|(d+*>-2_ (6 



n k. 



sing 

Fisher and Fisher [|o|] continued to argue that there are two alternatives. Either k ~ 
K s ing, implying zv — 1; or the full compressibility k is constant, implying z = d. The 
former is consistent with the Gaussian values v = | , z = 2 found by Uzunov p4| ] for 
the pure case in d < 2. The latter is believed to apply to repulsively interacting bosons 
in a random media. These remarkable simple arguments thus predict the exact value 
z = d for the dynamic exponent in this case. 

For later reference, let us consider the charged case and calculate the conductivity 
a. The only relevant term for this purpose is the first one in (^8) with Vy> replaced by 
Vip — eA. We allow the superfluid mass density to vary in space and time. The term 
in the action quadratic in A then becomes in the Fourier representation 



/ A(~fc ,-k)p s (fc ,k)A(fc ,k). (6.16) 

Jk ,U 



The electromagnetic current, 



^■ k » = <A(-t-k) (<U7) 
obtained from this action can be written as 

j(fc ,k) -a(fc ,k)E(fc ,k) (6.18) 

with the conductivity 

CT (fc) = ^^) (6 .i9) 
k 

essentially given by the superfluid mass density divided by ko, where it should be 
remembered that the mass m is effectively set to 1 here. 

The above hyperscaling arguments have been extended by Fisher, Grinstein, and 
Girvin | f77| ] to include the l/|x|-Coulomb potential. The quadratic terms in the effec- 
tive theory (3.73) may be cast in the general form 

4 2 ff } = ^(p S k 2 -^^Kfc)| 2 , (6.20) 

where e is the renormalized charge. From ( 3.73| ) we find that to lowest order: 

e 2 = 2 d -V d - 1 )/ 2 r[i(d-l)]eg. (6.21) 

The re norma lized charge is connected to the (0 0)-component of the full polarization 
tensor ( 3.19 ) via 

e 2 =lim ' ' ■ (6.22) 
|k|-o IT oo (0,k) 
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Table 6.1: The upper critical space dimension d c of a nonrelativistic (NR) and a rela- 
tivistic (R) quantum theory with a \4>\ 2k interaction term. 



k 


dc(NR) 


4(R) 


2 


2 


3 


3 


1 


2 


oo 





1 



A simple hyperscaling argument like the ones given above shows that near the transi- 
tion, the renormalized charge scales as 



e 



2 



^- Z . (6.23) 



They then argued that in the presence of random impurities this charge is expected to 
be finite at the transition so that z = 1. This again is an exact results which replaces 
the value z = d of the neutral system. 

We have seen that d c = 2 is the upper critical dimension of the nonrelativistic 
1 0| 4 -theory. Dimensional analysis shows that for an interaction term of the form 



the upper critical dimension is 



Ci = -goW (6.24) 



(6.25) 



k-1 



The two important physical cases are d c = 2, k = 2 and d c = 1, k = 3, while d c — > 
when k — > oo. For space dimensions d > 2 only the quadratic term, \(f>\ 2 , is relevant 
so that here the critical behavior is well described by the Gaussian theory. 

In the corresponding relativistic theory, the scaling dimensions of t and x are, of 
course, equal [t] = [x] = — 1 and [(f)] = \{d— 1). This leads to different upper critical 
(space) dimensions, viz., 

4 = ^ = ^ + 1, (6.26) 



instead of ( 6.25 ). The two important physical cases are here d c — 3, k — 2 and d c = 2, 
k = 3, while d c — > 1 when k — > oo. On comparison with the nonrelativistic results, we 
see that the nonrelativistic theory has an upper critical space dimension which is one 



lower than that of the corresponding relativistic theory (see Table 6.1 ). Heuristically, 
this can be understood by noting that in a nonrelativistic context the time dimension 
counts double in that it has a scaling dimension twice that of a space dimension [see 



Eq. (6.6)], thereby increasing the effective spacetime dimensionality by one. 

From this analysis it follows that for a given number of space dimensions the crit- 
ical properties of a nonrelativistic theory are unrelated to those of the corresponding 
relativistic extension. 
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In closing this section we recall that in a one-dimensional relativistic theory — 
corresponding to the lowest upper critical dimension (d c = 1) — a continuous sym- 
metry cannot be spontaneously broken. However, the theory can nevertheless have a 
phase transition of the Kosterlitz-Thouless type. Given the connection between the 
relativistic and nonrelativistic theories discussed above, it seems interesting to study 
the nonrelativistic theory at zero space dimension (d = 0) to see if a similar rich phe- 
nomenon as in the lower critical dimension of the relativistic theory occurs here. This 
may be of relevance to so-called quantum dots. 



6.3 Quantum Hall Liquid 

In this section we shall argue that the effective theory of a quantum Hall liquid can be 
used to describe its liquid-to-insulator transition as the applied magnetic field changes, 
and study its critical properties. 

Experimentally, if the external field is changed so that the filling factor vh moves 
away from an odd-denominator value, the system eventually becomes critical and un- 
dergoes a transition to an insulating phase. Elsewhere [p7[], we have argued that this 
feature is encoded in the CSGL theory. In the spirit of Landau, we took a phenomeno- 
logical approach towards this field-induced phase transition. And assumed that when 
the applied magnetic field H is close to the upper critical field H+ H at which the 
quantum Hall liquid with filling factor is destroyed, the chemical potential of the 
composite particles depends linearly on H, i.e., /i oc ei?+ — eH. This state can 
of course also be destroyed by lowering the applied field. If the system is near the 
lower critical field H~ , we assumed that the chemical potential is instead given by 
Ho oc eH — eH~ H . This is the basic postulate of our approach. 



We modify the CSGL Lagrangian (5.2) so that it only includes the fluctuating part 



of the statistical gauge field. That is, we ignore the classical part of a which yields a 
magnetic field that precisely cancels the externally applied field. We can again trans- 
form the mass m of the nonrelativistic |0| 4 -theory away. In addition to the engineering 



dimensions (6.6), we have for the Chern-Simons field 

[ea,] = l, [ea }=2, [9} = 0. (6.27) 

In two space dimensions, the coupling constant Ao was seen t o be dimensionless, im- 
plying that the |0| 4 -term is a marginal operator. From (6.27) it follows that also the 



Chern-Simons term is a marginal operator. Hence, the CSGL theory contains — apart 
from a random and a Coulomb term — precisely those terms relevant to the description 
of the liquid-to-insulator transition in a quantized Hall system. 

It is well-known [jO, ^3], that the coefficient of the Chern-Simons term is not renor- 
malized by quantum corrections. 

To first order in a loop expansion, the theory in two space dimensions was known 
to have an IR fixed point determined by the zero of the beta function [|8~3|l, 

(3(X) = - (iX 2 - 1 ) . (6.28) 



7T 
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Figure 6.5: Schematic representation of the beta function (6.28 ). 



The calculation of /3(A) has been extended to fourth order in the loop expansion [|l7j]. 
This study revealed that the one-loop result ( 6.28 ) is unaffected by these higher-order 
loops. Presumably, this remains true to all orders in perturbation theory, implying 
that — just as in the neut ral sy stem which corresponds to taking the limit 8 — > oo — the 
one-loop beta function (5.28) is exact. 



It is schematically represented in Fig. 5.5, and is seen to yield a nontrivial IR fixed 
point A* 2 = 1/40 2 determined by the filling factor. More precisely, the strength of the 
repulsive coupling at the fixed point A* = n(2l + 1) increases with the number 21 + 1 
of flux quanta bound to the electron. The presence of the fixed point shows that the 
CSGL theory undergoes a 2nd-order phase transition when the chemical potential of 
the composite particles tends to a critical value. As in the neutral case, it can be shown 
that the boson self-energy S also vanishes at every loop order in the charged theory, 
and that the self-coupling parameter A is the only object that renormalizes. The 2nd- 
order phase transition described by the nontrivial IR fixed point has consequently again 
Gaussian exponents v = |, z = 2 [17 1. It should be noted that only the location of 
the fixed point depends on 8, the critical exponents — which in contrast to the strength 
of the coupling at the fixed point are independent of the regularization and renormal- 
ization scheme — are universal and independent of the filling factor. This "triviality" is 
in accord with the experimentally observed universality of the FQHE. A dependence 
of the critical exponents on 8 could from the theoretical point of view hardly be made 
compatible with the hierarchy construction [ ]84| ] which implies a cascade of phase tran- 
sitions. From this viewpoint the present results are most satisfying: the CSGL theory is 
shown to encode a new type of field-induced 2nd-order quantum phase transition that 
is simple enough not to obscure the observed universality of the FQHE. 

We stress again that in order to arrive at a realistic description of the FQHE, the 
CSGL theory has to be extended to include a l/|x|-Coulomb potential and impuri- 
ties. Both will change the critical behavior we found in this section. In particular, the 
Coulomb potential will change the Gaussian value z = 2 into z = 1. 
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In Sec. 6.2, we saw that in the absence of impurities repulsively interacting bosons 
undergo a 2nd-order quantum phase transition as function of the chemical potential. 
As was pointed out there, this universality class is of relevance to various condensed- 
matter systems. However, in most of the systems mentioned there, as well in 4 He in 
porous media, impurities play an essential if not decisive role. For example, the two- 
dimensional superconductor-to-insulator transition investigated by Hebard and Paala- 
nen [ |63| ] is driven by impurities. This means that, e.g., the correlation length £ diverges 
as | A c — A | ~ v when the disorder strength A characterizing the randomness approaches 
the critical value A c . Hence, a realistic description of the critical behavior of these sys- 
tems should include impurities. 



To include these, we proceed as before and add the random term (3.75) to the 
nonrelativistic |0| 4 -theory (3.1). The random field ^( x ) has the Gaussian distribution 



(3.76). We shall study the theory in the symmetrical state where the bare chemical po- 
tential is negative and the global U(l) symmetry unbroken. We therefore set /i = — r 
again, with tq > 0. We leave the number of space dimensions d unspecified for the 
moment. As we remarked before, since tp(x) depends only on the d spatial dimensions, 
the impurities it describes should be considered as grains randomly distributed in space. 
When — as is required for the study of quantum critical phenomena — time is included, 
the static grains trace out straight worldlines. That is to say, these impurities are line- 
like in the quantum theory. It has been shown by Dorogovtsev [ ]86| ] that the critical 
properties of systems with extended defects must be studied in a double e-expansion, 
otherwise no IR fixed point is found. The method differs from the usual e-expansion, 
in that it also includes an expansion in the defect dimensionality ed- To carry out this 
program in the present context, where the defect dimensionality is determined by the 
dimensionality of time, the theory has to be formulated in ed time dimensions. The 
case of interest is ed = 1, while in the opposite limit, ed — > 0, the random nonrelativis- 
tic |</>| 4 -theory reduces to the classical spin model with random (pointlike) impurities. 
Hence, ed is a parameter with which quantum fluctuations can be suppressed. An ex- 
pansion in ed is a way to perturbatively include the effect of quantum fluctuations on 
the critical behavior. Ultimately, we will be interested in the case ed = 1. 

To calculate the quantum critical properties of the random theory, which have first 
been studied in |]87|], we will not employ the replica method [|8[, but instead follow 
Lubensky [|8f^]. In this approach, the averaging over impurities is carried out for each 
Feynman diagram separately. The upshot is that only those diagrams are to be included 
which remain connected when Ao, the parameter characterizing the Gaussian distribu- 
tion of the impurities, is set to zero [p0|]. To ob tain t he relevant Feynman rules of the 
random theory we average the interaction term (3.75) over the distribution (3.76): 



Dtp P( , 0) exp 



exp 



ii 2 *«A J d^dWM^i^x)! 2 !^',*)!' 



(6.29) 



The randomness is seen to result in a quartic interaction term which is nonlocal in time. 
The factor i £d appearing in ( |6l9| ) arises from the presence of ed time dimensions, each 
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of which is accompanied by a factor of i. The Feynman rules of the random theory are 
now easily obtained 



_j-«dgi(wi+W2H hw £d )?7 



LUl + LU-2 H hW Ed 

-4i £d A 



k 2 — rg + i?/ 



= ^ ed (2^) ed ( 5 ed (c i ; 1 +cJ 2 + •••+cc )ed )Ao, 



(6.30) 



where we note that the Lagrangian in e d time dimensions involves instead of just one 
time derivative, a sum of e d derivatives: d t d tl + d t2 + ■ ■ ■ + d tcd ■ 

Following Weichman and Kim [pTp, we evaluate the integrals over loop energies 
assuming that all energies are either positive or negative. This allows us to employ 
Schwinger's properti me re presentation of propagators [^2B, which is based on the in- 
tegral representation ( 3.44 ) of the gamma function. The energy integrals we encounter 
to the one-loop order can be carried out with the help of the equations 



1 



(27r) £d ui\ + u>2 

r(i - e d ' 



(2tt)« 



sgn(x)|cc| ed 

p i(aJi+0J 2 H 



d - x ± ir\ 
1 ^ e ±isgn(a)7re d _|_ ^ 



(27T) £ d UJi + UJ 2 H hi 

■my*- 1 



IXT) 



(27r)^r(e d 



sin(±7re d ) 



sgn(a;) 



sin( 5 7red) 



(6.31) 



(6.32) 



The prime on the integrals is to remind the reader that the energy integrals are taken 
over only two domains with either all energies positive or negative. Th e ene rgy in- 
tegrals have been carried out by using again the integral representation ( 3.44 ) of the 
gamma function. In doing so, the integrals are regularized and — as is always the case 
with analytic regularizations — irrelevant diver genc es suppressed. 

By differentiation with respect to x, Eq. ( |5.3l| ) can, for example, be employed to 
calculate integrals involving integrands of the form l/(k>i +ci>2 + - ■ ■+o->e d —x + ir)) 2 . Is 
is easily checked that in the limit e d — * 1, wh ere th e ener gy integral can be performed 
with help of contour integration, Eqs. ( 6.31 ) and ( J6.32 ) reproduce the right results. 
When considering the limit of zero time dimensions (e d — > 0), it should be borne in 
mind that the energy integrals were taken over two separate domains with all energies 
either positive or negative. Each of these domains is contracted to a single point in 
the limit e d — * 0, so that one obtains a result which is twice that obtained by simply 
purging any reference to the time dimensions. The integral ( 5.32) contains an additional 
convergence factor exp(iw^) for each of the e d energy integrals. This factor, which — 
as we remarked before — is typical for nonrelativistic quantum theories [pT[], is to be 
included in self-energy diagrams containing only one ^-propagator. 
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Before studying the random theory, let us briefly return to the repulsively inter- 
acting bosons in the absence of impurities. In this case, there is no need for an ej- 
expansion and the formalism outlined above should yield results for arbitrary time 
dimensions < ea < 1, interpolating between the classical and quantum limit. After 



the energy integrals have been performed with the help of Eqs. (5.31) and (5.32), the 



standard technique of integrating out a momentum shell can be applied to obtain the 
renormalization-group equations. For the correlation-length exponent v we obtain in 
this way Q 



1 

v = — 
2 



1 



2 (to + 4) - (to + 3)e d 



C0S 2 (57T£d) 



(6.33) 



Here, e = 4 — 2ed — d is the deviation of the effective spacetime dimensionality from 
4, where it should be noted that in (canonical) nonrelativistic theories, time dimen- 
sions have an engineering dimension twice that of space dimensions. (This property is 
brought out by the Gaussian value 2 = 2 for the dynamic exponent z.) For compar- 
ison we have extended the theory (3.1) to include to complex ^-fields instead of just 
one field. In the classical limit, Eq. ( 5.33| ) gives the well-known one-loop result for a 
classical spin model with 2m real components [29|, 



1 



e to + 1 
2 to + 4 



(6.34) 



while in the quantum limit it gives the result u — > |, as required. 

The exponent (5.33), and also the location of the fixed point, diverges when the 
number of time dimensions becomes ed — > (m + 4)/(m + 3). Since this value is 
always larger than one, the singularity is outside the physical domain < e<j < 1. This 
simple example illustrates the viability of the formalism developed here to generate 
results interpolating between the classical and quantum limit. 

We continue with the random theory. After the energy integrals have been carried 
out, it is again straightforward to derive the renormalization-group equations by inte- 
grating out a momentum shell A/6 < k < A, where A is a high-momentum cutoff and 
b = cxp(Z), with I infinitesimal. Defining the dimensionless variables 



where 



A=-^AA- e : 
(27r) £ i 



A = K d AA 



d-4. 



= rA~ 



(4^) d / 2 r(f) 

is the area of a unit sphere in d spatial dimensions divided by (27r) d , we find [E 



(6.35) 



(6.36) 



dA 
dZ 
dA 
"dT 
df 
d7 



eA - 8 [r(l - e d ) + (to + 3)r(2 - e d )] cos(i7re d )A 2 + 6AA 



— = (e + 2e d )A + 4A 2 - 16(TO + l)r(2-e d )cos(i7red)AA 



2f + 47T 



m + 1 cos 2 (i7red) 
r(e d ) sin(i7re d ) 



A- A. 



(6.37) 
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These results are to be trusted only for small values of e d . For illustrative purposes we 
have, however, kept the full e d dependence. The set of equations yields the fixed point 



A* 



A* = 



1 



6e d 



16cos(±7re d )r(l - e d ) 2m(l - e d ) - 1 

l m(l-e d )( 2e d - e) + 2e d (4 - 3e d ) + e(2 - e d ) 
4 



(6.38) 



and the critical exponent 
1 

V =2 + 



2m(l - e d ) - 1 
e + 2e d m + 1 (6e d + e) [e d + cos(7re d )] 



16 



16 



2m(l - e d ) - 1 



(6.39) 



The dynamic exponent is given by z = 2 + A*. When the equations are expanded to 
first order in e d , we recover the IR fixed point found by Weichman and Kim [ ^l| ] using 
an high-energy cutoff: 



A* = 



16 2m 

with the critical exponent 



1 

V= 2 



1 e + 6e d . ^ _ 1 (2-m)e + 2(m + 4)e d 
1 ' ~~ 4 2m - 1 



1 + 



1 3me + (5m + 2)2e d 



2m- 1 



(6.40) 



(6.41) 



The value of the critical exponent (5.41 ) should be compared with that of the classi- 
cal spin model with 2m components in the presence of random impurities of dimension 

1 3me + (5m + 2)e d 



1 



2m- 1 



(6.42) 



Taking into account that in a nonrelativistic quantum theory, time dimensions count 
double as compared to space dimensions, we see that both results are equivalent. As to 
the dynamic exponent, notice that the perturbative result z = 2 + A*, with A* given 
by ( |5.40 ), is far away from the exact value z = d for e d = 1 [p0[|. 



6.5 Experiments 

6.5.1 Superconductor- To-Insulator Transition 

The first experiments we wish to discuss are those performed by Hebard and Paalanen 
on superconducting films in the presence of random impurities [ ^3[ p5j]. It has been 
predicted by Fisher that with increasing applied magnetic field such systems un- 
dergo a zero-temperature transition into an insulating state. (For a critical review of the 
experimental data available in 1993, see Ref. [0].) 

Let us restrict ourselves for the moment to the TA-plane of the phase diagram by 
setting the applied magnetic field H to zero. For given disorder strength A, the system 
then undergoes a Kosterlitz-Thouless transition induced by the unbinding of magnetic 
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vortex pairs at a temperature Tkt well below the bulk transition temperature (see Sec. 



4.3 ). The Kosterlitz-Thouless temperature is gradually suppressed to zero when the 
disorder strength approaches criticality A — > A c . The transition temperature scales 
with the con-elation length £ — |A C - A|~^ as Tkt ~ £,~ z - 

In the i/A-plane, i.e., at T = 0, the situation is as follows. For given disorder 
strength, there is now at some critical value H c of the applied magnetic field a phase 
transition from a superconducting state of pinned vortices and condensed Cooper pairs 
to an insulating state of pinned Cooper pairs and condensed vortices. The condensa- 
tion of vortices disorder the ordered state as happens in classical, finite temperature 
superfluid- and superconductor-to-normal phase transitions [ f34| ] . When the disorder 
strength approaches criticality again, H c is gradually suppressed to zero. The critical 
field scales with £ as H c ~ $o/£ 2 - in fact, this expresses a more fundamental result, 
namely that the scaling dimension c?a of A is one, 

d A = 1, (6.43) 

so that |A| — From this it in turn follows that E ~ t,^ 1 ^ 1 — £~ (z+1) , and that 
the scaling dimension d,Ao °f A) ls z, 

d A „ = z, (6.44) 

so that Aq ~ £ t _1 ~ £~ z . Together, the scaling results for Tkt and H c imply that 

H C ~T^. (6.45) 

This relation, linking the critical field of the zero-temperature transition to the Kosterlitz- 
Thouless temperature, provides a direct way to measure the dynamic exponent z at the 
H = 0, T = transition. This has been first done by Hebard and Paalanen [^3| |85| ]. 
Their experimental determination of Tkt and H c for five different films with varying 



amounts of impurities confirmed the relation (6.45) with 2/z = 2.04 ± 0.09. The zero- 
temperature critical fields were obtained by plotting dp xx / dT | h versus H at the lowest 
accessible temperature and interpolating to the the field where the slope is zero. The 
resulting value z — 0.98 ± .04 is in accordance with Fisher's prediction [p3|, z = 1, 
for a random system with a l/|x| -Coulomb potential. 



Hebard and Paalanen [ 63 1 also investigated the field-induced zero-temperature tran- 
sition. The control parameter is here S oc H — H c . When plotted as function of 
\H - H C \IT X I VHZH they saw their resistivity data collapsing onto two branches; an 
upper branch tending to infinity for the insulating state, and a lower branch bending 
down for the superconducting state. The unknown product vhzh is experimentally 
determined by looking for which value the best scaling behavior is obtained. Further 
experiments carried out by Yazdani and Kapitulnik [ j?4] ] also determined the product 
vh{zh + 1) (see below). The two independent measurements together fix the crit- 
ical exponents vh and zh separately. From their best data, Yazdani and Kapitulnik 
extracted the values [Q] 

z ff = 1.0 ± 0.1, v H = 1-36 ± 0.05. (6.46) 
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6.5.2 Quantum-Hall Systems 

We continue to discuss the field-induced quantum phase transitions in quantum Hall 
systems. Since an excellent discussion recently appeared in the literature [Q], we shall 
be brief, referring the reader to that review for a more thorough discussion and addi- 
tional references. 

One can image transitions from one Hall liquid to another Hall liquid with a dif- 
ferent (integer or fractional) filling factor, or to the insulating state. Experiments seem 
to suggest that all the quantum-Hall transitions are in the same universality class. The 
transitions are probed by measuring the resistivities p xx and p x v . From the dependence 
of the conductivity a on the superfiuid mass density, Eq. ( 6.19] ), and the scaling relation 



( |6.14[ ), it follows that it scales as [|95|] 

a~C {d ~ 2) - (6.47) 

In other words, the scaling dimension of the conductivity and therefore that of the 
resistivity is zero in two space dimensions. On account of the general finite-size scaling 



form (5.3), we then have in the limit |k| — > 0: 

p xx/y (k , H, T) = g xx/y (k /T, \Sr/T), (6.48) 

where the distance to the zero-temperature critical point is measured by S oc H — 
H^ H ~ t x I vz . This scaling of the width of the transition regime with temperature has 
been corroborated by DC or k = experiments on various transitions between integer 
quantum-Hall states which were all found to yield the value 1/vz — 0.42 ± 0.04 [p6|. 

A second measurement of the critical exponents involves the applied electric field. 
As we have seen above, it scales as E ~ so that for the DC resistivities we 

now obtain the scaling form: 

Pxx/y (H, T, E) = g xx/y (\Sr/T, \8\ V ^/E), (6.49) 

The scaling \5\ ~ tfV^+i) has again been corroborated by experiment which yielded 
the value v(z + 1) w 4.6 Together with the previous result obtained from the 
temperature scaling this gives 

z«l, i/w2.3. (6.50) 

The value of the dynamic exponent strongly suggests that it is a result of the presence 
of the l/|x|-Coulomb potential. The correlation length exponent v is large. 

6.5.3 Id Electron Systems 

Recently, sil icon MOS FET's at extremely low electron number densities has been stud- 
ied [ p8[ p9 , 100| , 101]. Earlier experiments at higher densities seemed to confirm the 



general believe, based on the work by Abrahams et al. [^5j|, that such two-dimensional 
electron systems do not undergo a quantum phase transition. In that influential paper, 
it was demonstrated that even weak disorder is sufficient to localize the electrons at the 
absolute zero of temperature thus excluding conducting behavior. Electron-electron 
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interactions were however not included. As we saw in Sec. 3.2, at low densities, the 
l/|x| -Coulomb interaction becomes important and the analysis of Abrahams etal. [ |9^ ] 
no longer applies. 

The recent experiments have revealed a zero-temperature conductor-to-insulator 
transition triggered by a change in the charge carrier density n. That is, the distance to 
the critical point is in these systems measured by S oc n — n c . Like in the quantum- 
Hall systems, these transitions are probed by measuring the resistivity. It scales with 



temperature near the transition according to the scaling form (5.49) with H set to zero. 
For n < n c , where the Coulomb interaction is dominant and fluctuations in the charge 
carrier density are suppressed, the electron system is insulating. On increasing the 
density, these fluctuations intensify and at the critical value n c , the system reverts to a 
conducting phase. By plotting their conductivity data as function of T / \ 8 \ vz with vz = 



1.6 ± 0.1, Popovfc, Fowler, and Washburn [101] saw it collapse onto two branches; 
the upper branch for the conducting side of the transition, and the lower one for the 
insulating side. A similar collapse with a slightly different value X/vz = 0.83 ± 0.08 



was found in Ref. [99|, where also the collapse of the data when plotted as function of 
\ s y E i/v{z+i) was obtained. The best collapse resulted for l/(z + l)v = 0.37 ± 0.01, 
leading to 

« = 0.8 ±0.1, i/=1.5±0.1. (6.51) 

The value for the dynamic exponent is close to the expected value z = 1 for a charged 
system with a l/|x|-Coulomb interaction, while that of v is surprisingly close to the 



value (6.46) found for the superconductor-to-insulator transition. 

A further experimental result for these two-dimensional electron systems worth 
mentioning is the suppression of the conducting phase by an applied magnetic field 



found by Simonian, Kravchenko, and Sarachik [ 100|. They applied the field parallel 
to the plane of the electrons instead of perpendicular as is done in quantum-Hall mea- 
surements. In this way, the field presumably couples only to the spin of the electrons 
and the complications arising from orbital effects do not arise. At a fixed temperature, 
a rapid initial raise in the resistivity was found with increasing field. Above a value of 
about 20 kOe, the resistivity saturates. It was pointed out that both the behavior in a 
magnetic filed, as well as in zero field strongly resembles that near the superconductor- 
to-insulator transition discussed above, suggesting that the conducting phase might in 
fact be superconducting. 



6.5.4 Conclusions 

We have seen that general scaling arguments can be employed to understand the scal- 
ing behavior observed in various quantum phase transitions. Most of the experiments 
seem to confirm the expected value z = 1 for a random system with a l/|x|-Coulomb 
interaction. The number of different universality classes present is yet not known. Even 
if the conductor-to-insulator transition observed in silicon MOSFET's at low electron 
number densities turns out to be in the same universality class as the superconductor- 
to-insulator transition, there are still the field-induced transitions in quantum-Hall sys- 
tems, which have a larger correlation-length exponent. 

The paradigm provided by a repulsively interacting Bose gas, seems to be a good 
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starting point to describe the various systems. However, high-precision estimates cal- 
culated from this theory with impurities and a l/|x|-Coulomb interaction included are 
presently lacking. 
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